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ABSTRACT 


The external nucleate росі воля heat-transfer 
coefficient of both smooth and enhanced horizontal tubes 
in R-1i14-o1l mixtures (O to 10 mass percent oil) were 
measured for heat fluxes from 1 to 91 kW/m? at a saturation 
temperature о? 2.2 С. The enhanced tubes tested were a 
Wieland GEWA-T finned tube containing 1.02 fins/mm, a 
Hitachi  Thermoexcel-E and a Hitachi  Thermoexcel-HE tube. 
The Thermoexcel-E and -HE tubes with their re-entrant 
cavity designs were found to improve the heat-transfer 
coefficient over the smooth tube value at a constant heat 


flux by» a factor of approximately y in oil-free R iR 


` while the GEWA-T tube improved the coefficient by a factor 


of about 4. While all of the tubes showed a generally 
decreasing performance with the presence of oil. the GEWA-T 
tube resulted in the minimum reduction. This tube showed 


up to 20% and 35% reduction (compared to the oil-free case) 
at 3 percent and 10 percent oi respectively. The 


Thermoexcel-E and Thermoexcel-HE tubes showed performance 


reductions of up to 40* апа 60%, respectively, with W 
percent oil. At a practical heat flux of about 30 kW/m° 
with 3 percent o. the СЕМА-Т, Thermoexcel-E and 
Thermoexcel-HE tubes outperformed the smooth tube Бу 


factors of 4.8, 4.6 and 4.0, respectively. 
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INTRODUCTION 


ast + на = + 


m BACKGROUND 

In the design of refrigeration systems for use aboard 
Naval vessels, the desire is to make the System as 
efficient as possible. To achieve this higher efficiency, 
one must consider the use of a high pressure refrigerant, 
since this gives a higher energy transfer per unit volume 
ЕШ уарог (і.е.. Ey = Reg /v,) (Ref. 11], and results in less 
of a pressure drop along the vapor path for a given heat 
load. Го, which is а high pressure refrigerant, has an 


energy transfer per unit volume of vapor of 8570 kJzm* 


(Ref. 21 апа, since it is higher pressure, it requires 
heavier gauge components. E dul. on the other hand, is a 
low pressure refrigerant witha low energy transfer per 
unit volume of vapor of 2086 kJ/m? and requires lighter 


gauge material. 

Table 1 provides a short summary of some of the 
available refrigerants and their respective energy transfer 
per unit volume of vapor, their relative toxicity and their 
required system pressure for optimum operation. As can be 
seen from Table 1, R-12 and R-22 possess the highest values 
of energy transfer per unit volume of vapor, and therefore 
would provide systems of higher efficiency than the other 
refrigerants listed since pressure losses in the system 


would be less. However, since both R-12 and R-22 are 


11 


high-pressure refrigerants, their respective systems would 

be significantly heavier larger than those of the-other 

refrigerants. However, R-114 being a moderate-pressure 

refrigerant with an energy transfer per unit voulume of 

vapor of 3800 КЈУШУ, ӘС смеше allow for а relatively 

lighter component system than those of Н-12 апа 2 
TABLE 1 Ret. wed 


SUMMARY OF REFR1GERANTS 
Operational 


System 

Refrigerant ..К43/П! |. Тохїс [ыыы ИИ 
к=! 2086 moderate (Group 5D low 
Ниша 8570 low (Group 6) high 
него 12480 moderate (Group S) high 
к-113 1125 high (Group 4) low 

K-14 3800 moderate (Group 6) moderate 


Additionally, R-114 belongs to the refrigerant group (Group 


6 ) having the lowest toxicity [ReF. 21, апа it has also 
been proven fairly stable with temperature [ Ref. 23” 
Therefore, R-114 is а viable alternative to the higher 


pressure refrigerants. 


Ва REVIEW OF REFRIGERANT-OIL MIXTURE BOILING BEHAV GER 

The presence of ea in the evaporators Of 
refrigeration systems is a common occurrence. This is due, 
in~ part, to -the Use nek lubricants in the hermetically 
sealed refrigeration compressors. Over a period of time, a 
considerable amount of өлі is introduced inte the 
refrigerant, subsequently changing the physical 
properties of the refrigerant. Jensen and Jackman [Ref. 3) 
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reported that the density and the specific heat behave 
linearly іп the refrigerant-o1! mixtures, however, the 
viscosity and the surface tension do not. Henrici and 


Hesse [Ref. 4] experimentally determined in 1971 that the 


surface tension of mixtures first decreased for R-114-oil 
mixtures of up to 2.5 percent oil, and at higher oil 
concentrations, the surface tension increased. This non- 


linear behavior of the properties of the mixtures makes it 
difficult to explain the changes in the heat-transfer 
coefficient for these mixtures. 

The most evident example of the addition of oil to the 


refrigerant is the presence of foam in the mixture when 


boiling takes place. With the oil concentration above 1 
percent, there is a substantial amount of foam generated 
from the nucleate boiling surface. This foam is generated 


because the refrigerant in the refrigerant-oil mixture is 
more volatile than the oil and therefore can vaporize, 
creating a gas bubble surrounded Бу ап oil-rich layer 


(Кет. 5]. Figure 1.1 shows the ideal model of the growth of 


a bubble in a refrigerant-oil mixture. These generated 
bubbles are at lower density than that of the surrounding 
refrigerant-oil mixture. The bubbles ascend to the liquid- 


vapor interface and collect at the interface and produce a 
foam layer. As reported by Henrici and Hesse (Ref. 4], the 


foaming action is most pronounced іп refrigerant-oil 
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mixtures of БӘРЕП —percentw"o11, and provide а the 
presence of oil results in’ a subsequent decrease in the 
heat-transfer performance of the boiling tube. 


Numerous researchers have tried to explain the reasons 


for the decrease in the heat-transfer coefficient in 
refrigerant-oil mixtures of increasing oi! content. Thome 
[ Ref. 6], after an extensive literature search, found that, 


in 1956, Van Wijk et al. explained the effect as being the 
result of the evaporation of the more-volatile component 
leaving an oil-rich layer which has a higher resultant 
local boiling point. This increased local boiling point 
requires an increase in the amount of superheat necessary 
to continue the vaporization and bubble growth. The result 
is a decrease in the heat-transfer coefficient. 


Stephan and Preusser (Ref. 7] determined that the work 


required to form bubbles in the refrigerant-oil mixtures 
was greater than in an equivalent amount of pure 
refrigerant. Therefore, the heat-transfer coefficient was 


lower for the refrigerant-oil mixture than for the pure 
refrigerant. Chongrungreong and Sauer (Ref. 8] reported 


that the rate of heat diffusion, which is governed by the 


thermal properties of the oil-rich layer, limits the bubble 
growth and that the effects of Surface tension are 
negligible. Thome [Ref. 61] agreed with the above factors 


and concluded that the viscosity variations are important 
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in explaining the decrease in the heat-transfer 
coefficient for some refrigerant-oil mixtures. 
Although researchers have not decided on any one set 


Of factors that determine the reasons for the decrease in 


the  heat-transfer coefficient as oil content increases in 
refrigerant-oil mixtures, all agreed that the physical and 
thermal properties are important factors in explaining the 


heat-transfer behavior of refrigerant-oil mixtures. 


eu REVIEW OF BOILING PERFORMANCE FROM SMOOTH TUBES 

пе Ue. oS; Navy attempts to reduce the size of Tes 
Shipboard refrigeration systems by the use of alternate 
refrigerants (such as R-114). At the same time, the Navy 
desires to increase the performance of these systems by the 
use of enhanced surfaces. However, as discussed by 
Wanniarachchi et al. (Ref. 9], only a limited amount of 
data is available regarding the boiling of refrigerants in 


the presence of oil. 


Stephan (Ref. 101 conducted experiments using a 
smooth horizontal plate as a boiling surface "in КЕЕ 
and R-22 refrigerants and reported up to 50 percent 
reductions in the boiling heat-transfer coefficient 
with a refrigerant-oil mixture of 9 percent oil (by 
mass), while 50 percent ou resulted in a 90 


percent reduction in the boiling heat-transfer coefficient. 
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In 19715; Henrici and Hesse [Ref . 4] sonauc ted 
experiments involving the boiling of R-114 in the presence 
2, oil from a smooth copper tube and reported a maximum 
decrease of 20 percent decrease in the boiling heat- 
transfer coefficient with 1 percent oil] by mass. b heed, 
Stephan and Mitrovic (Ref. 11] reported that for the GEWA-T 
tube (manufactured Бу the Wieland Company) Phe 1184-01 1 
mixtures, the boiling heat-transfer coefficient of the 
enhanced surface was significantly altered by the presence 
ere the oil. ШӘ Т wet lly (Ref. 5] reported à drop in 
the  heat-transfer coefficient of O to 35 percent (depending 
on the heat flux) for a R-114-oil mixture of 10 percent oil 


for a smooth copper tube. 


p. REVIEW OF BOILING PERFORMANCE FROM ENHANCED SURFACES 
Arai et al. (Ref. 121 tested a 200-ton. 5-12 
centrifugal water chiller. The chiller which used ап 


enhanced-surface tube, commercially known as Thermoexcel-Em 


made by Hitachi, was 28 percent shorter in length and had 
an approximate improvement in the overall heat-transfer 
coefficient of SO percent. Additional research by Yilmaz 


and Westwater CRef. 1951, Marto and Lepere (Ref. 141, and 
Carnavos (Ref. 151 on commercially available enhanced 
surface tubes, in various refrigerants other than R-114, 
indicated that a porous-coated surface exhibited the best 


НӨ пе heat-transfer coefficient ina pure refrigerant. 


I 


Comparisons by Reilly [Re tT] indicated that Sin R 2000 
Union Carbide's "High Flux" tube with а регоц5 coaten 
surface resulted in an improvement of at least a factor of 
seven over a smooth copper tube at the same operational 
parameters. Reilly's findings complemented the research of 
Yilmaz and Westwater (Ref. 112 Additionally, Reilly 
reported that for a porous-coated surface in an R-114-oil 
mixtures of 1 percent oil, the oil was seen to degrade the 
nucleate pool-boiling heat-transfer performance by about 
2O percent. 

16 Surfaces with Reentrant Cavities 

Reentrant cavities have been investigated for a 
considerable period of time. Griffith and Wallis (Ref. 151 
determined that the reentrant cavity geometry is important 
in the generation of nucleate boiling sites in two ways; 
Gi) the diameter of the reentrant cavity mouth determines 
the amount of superheat required CO initiate boiling, 
and (020) the cavity shape determines the stability once 


the boiling has begun. Webb (CRef. 16] stated that the 


key to the high performance of the reentrant grooved 
structures can be attributed to three factors: ( 
reentrant cavity within a “~€rirtical size range, (2) 


interconnected cavities, and (3) nucleation sites of a re- 
entrant shape. If the cavities are interconnected, 


adjacent cavities can activate each other. Webb (Ref. 16] 


15 


also stated that reentrant cavities provide a stable vapor 
trap, which can remain active at low values of superheat. 

The | GEWA-T (1.02 fins/mm)  enhanced-surface tube, 
manufactured by Wieland-Werke AG, is a copper tube with 
spiral fins of a T-shaped profile. It has a reentrant-type 
cavity formed by cold working an integral finned tube. The 
two adjacent T-fins form a spiral cavity with an opening 
of O.25-mm (see Figure 1.2). 

The Thermoexcel-E and Thermoexcel-HE enhanced-surface 
tubes are made by the Hitachi Company. Like the GEWA-T 
tube, the Thermoexcel-E (Figure 1.3) and Тћегтоехсе!-НЕ 


(Figure 1.4) tubes are also formed by cold working integral 


fins; however, these tubes have low fins which have a 
small Space cut out at the fin tips, giving the appearance 
ОҒ a sawtooth. These sawtooth fins are bent parallel to 
the tube axis, toran horizontal position, forming tunnels 
with evenly spaced surface pores. This cold working 


procedure results in a high area density of reentrant 
nucleation sites. The Thermoexcel-HE tube is a variation 
of the Thermoexcel-E tube, designed for improved 
performance in low-heat-flux regions. 

Experimental data showing the effect o f oil 
concentrations on the heat-transfer coefficient of the 
GEWA-T, the Thermoexcel-E and Thermoexcel-HE tubes in R-114 


are lacking. In an effort to gain more understanding of 
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(a) Schematic of GEWA-T Surface Cross Section 
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Schematic, Cut-away Cross Section and 


Micrograph of the GEWA-T Surface. 


20 


`a. араз Me 92 ЖУ "жш eh ТЕ АМУ 
АЛЫК А у 


T E oM m 4 ч бү, š - тоо б 
vB AOI T E, 2222 A ANSA yr, T. : мж 
72: ARES cu еу С: EGA хе : Ld 7 со а чи ја 
Ч р? ссх $. x 7%) 2. 2 XA йл n 
p on OS ENS Pop X € е 55$ Ken oe БЕ ^ AUS ж 
еН УА ЗУТ ROE, D 
4 T 7 7 b 


ere пет рај 


P >. KAS ^ ó 
рубок туит “о 
АЯ 


MANN = ~ ; 
- Dose x Ae Же ы 
орыг дэгээ ЭС 
~, МЕ an й өв 
Son Fag че C VAS 22. 
т ~. о хэв sy, Р 

* б л ~ ка ЭГ: E, 2 = 3 2 у Я aas 
ж A Im e Вы et: АЖ va m мы 
> хээд BSS к НЕВЕРНЫЙ m E M ова ча m 49 
2 YQ Мл ме 7, me м € 2“ š ^ 5 Е VEM s 
ои ИД Неа NAE аи SRT IY. ХХ VES са ims 
MK CO S iq a Ne RRR ок и 
OPES Й ERR OPO А мА A Ne g SOS AE ESS a он 
У ФС, ЗЫ BF NSIT De. онла a ER АА 
уч УУ A CREME Nt SiS Ср у Ce m 


~ 
> 


Уз 
м . 


PARA. 
З ` 
š f ох Ў 

А" 

Rowe ados. у си ae dsl? 


mt ND У n a M 


TR a. 
б 
GS < 
xo 
У 


vo es 


ST ITI cron риса ние 
udi 612249 2 НЕ ету 
Кесем еқ СА М Ж учете, 





Бійемге 1.3 Surface Photograph and Enlarged Outer 


Surface of the Thermoexcel-E Surface. 
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1.4 Surface Photograph and Enlarged Outer 


Surface of the Thermoexcel-HE Surface. 


the pool-boiling performance of these tubes, the David W. 
Taylor Naval Ship Research and Development Center sponsored 
this investigation. The smooth copper tube and the GEWA-T 
tube were provided by the Wieland Company. The Thermoexcel - 
E and Thermoexcel-HE tubes were provided by the Hitachi 
Company. 
EL BEESISCOBJECTIVE 
The objectives of this thesis are: 
dius take data on a smooth tube with the use of 
auxiliary heaters together with a convection shield 
to enable the collection of more data with a minimum 
of operator effort, 
=. take boiling data ona smooth copper tube in R-114 
to provide a baseline for follow-on comparison 
for enhanced tubes at oil concentrations of O, 1, 
2, З, 6 and 10 percent (by mass) аа boiling 
temperature of 2.2 C and a heat-flux range from 500 


wet to 91 kW/m*. 


© take data оп the following tubes for the above- 
-mentioned conditions: 


а. GEWA-T (1.02 fins/mm), 
b. Thermoexcel-E, and 


C. Thermoexcel-HE. 
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AS OVERALL APPARATUS 

The equipment used for this study consisted of nine 
basic components: (a) a Ругех- аа tee for the ров 
boiling of the R-114 liquid; (b) a Pyrex-glass tee for the 
condensing of the R-114 vapor; (c) an R-114 liqum 
reservoir; (d) à water-ethylene glycol mixture sump; (e) an 
R-12 refrigeration system with an external proportioning 
valve; (1) Gui reservoir; Cg) a vacuum pump; 
(h) instrumentation and data-acquisition system; and (i) a 
convection shield. Figure 2.1 shows schematically the 
arrangement of the various components. Complete details of 
the design, construction and operation of the apparatus 
are provided by Karasabun [Ref. 171, Reilly [Ref. 5] 
апа Wanniarchchi et al. (Ref. 1]. The original design of 
Karasabun was modified by the inclusion of the convection 


shield in the pool-boiling section as shown in Figures 2.2 


апа 2 53 

The operation of the apparatus consisted of the 
boiling of R-114 liquid in the lower glass tee (1), 
and condensing of the ва уарог in the upper 
glass tee соо The resultant  condensate was fed to 
the distribution tube within the boiling section by 
gravity. A water-ethylene glycol mixture from a 3O-gallon 
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Bi eune 2.3 Photograph of Convection Shield. 
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sump was used to cool the copper condensing coll кетер” 
condensing the R-114 vapor. This water-ethylene glycol 
mixture sump was maintained between -18 and -14 C by a 
one-half ton R-12 air-conditioning system. The mixture was 
pumped through the condensing coil and controlled manually 
with valve VC. 

Valve VC controls the amount of the water-ethylene 
glycol mixture to the condensing coil, thereby controlling 
the pressure within the test loop. For example, the opening 
of valve VC allows more of the water-ethylene glycol 
mixture to pass through the condensing coil, resulting in 
additional condensing of the R-114, and subsequently 
lowering the temperature and pressure within the entire 
system. 

As the subcooled R-114 condensate enters the boiling 
section of the apparatus through the distribution tube, it 
strikes the convection shield. This shield provides 
increased mixing of the subcooled R-114 condensate with the 
remaining R-114 Cor R-114/oil mixture) in the lower portion 
of the boiling section. The shield deflects the incoming 
condensate away from the boiling surface and allows for 
further mixing of the condensate апа the existing R-114 
liquid by means of the auxiliary cartridge heaters. This 
deflection of the incoming condensate allowed for a more 
even temperature distribution of the liquid R-114 and 


therefore decreased the previously reported temperature 
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gradients within the liquid region. This additional mixing 
Ore the liquid  R-114, coupled with the ЕУ to use 
auxiliary cartridge heaters in the liquid and allowed for 
the timely achievement of steady-state conditions without 
the previously observed thermal gradients in (һе liquid 
region of the tests done without the aid of the convection 
shield. It was also observed throughout the operation of 
the apparatus that the data obtained using the convection 
shield had a higher degree of repeatability than the data 


obtained without the shield. 


B. BOILING-TUBE CONSTRUCTION 

The tubes tested during the course of this study were 
the smooth copper tube, the GEWA-T (1.02 -fins/mm), the 
Thermoexcel-E and the Thermoexcel-HE tubes. Alil of these 
tubes had an active test и (1.е., heated) of length 
ої Eos 2-mm (8 in), with an outside diameter of 15.9 mm 
(5/8 in) and an inside diameter of 12.7 mm (1/2 in). The 
remaining 114.3-mm section on either’ side of the heated 
test section was smooth and unheated and did not nucleate 
under any heat flux or any oil concenteation condition. 
Figure 1.2 shows the schematic of the GEWA-T finned 
surface. Karasabun [Refs 171 described the analysis used 
in the data-reduction program, allowing for the two end 


surfaces to be treated as extended fins from the center 


section and subsequently accounting for the heat losses. 


The center section of each tube tested was heated by a 
1000-Ы 240-V stainless’ steel cartridge heater inserted in 
the center section of each tube. The cartridge heater was 
inserted in the boiling tube and surrounded by a copper 
sleeve with eight 1.3 mm by 1.3 mm thermocouple channels. 
Figure 2.4 shows the layout of the installed сиве тата 
boiling section. Figure 2.5 shows the details of the 
thermocouple layout. The thermocouple hot junctions were 
inserted into the copper sleeve and the copper sleeve was 
then tinned; and heat was maintained to keep the solder in 
a molten state as the sleeve was inserted and positioned in 
the center section of the boiling tube. Upon completion of 
the thermocouple installation, calibration oT the 
thermocouples was’achieved. Reilly [Ref. 5] determined in 
the course of his investigation that since the data- 
reduction program utilized the differences between the 
thermocouples in all computations, such as wall 
temperatures minus saturation temperatures, the corrections 
provided by the thermocouple calibration меге basically 
unnecessary since the calibration procedures were only 
necessary for items dependent on absolute temperature. 
Appendix A summarizes the procedure used by Karasabun 
(Ref. T and Reilly ( Ref. S] to calibrate the 


thermocouples. 
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(a) View of the boiling tube thermocouple locations as seen 
From the front of the experimental apparatus. 
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(b) Thermocouple sleeve unwrapped (at section А-А) to show 
the relative locations of the thermocouple channels 
(all dimensions in millimeters). 


LEFT-END VIEW RIGHT-END VIEW 





(c) End views of the boiling tube. 


Figure 2/5 Boiling Tube Thermocouples Channels. 
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С. CONVECTION SHIELD 


Initial tests on the smooth copper tube revealed an 
axial thermal gradient in the liquid over the length of the 
boiling tube. This thermal gradient increased with 


increasing oil concentration and with increasing heat flux, 
with the thermal еее тар! пр from 0.5 К to 1.5 K over 
the length of the boiling section. This thermal gradient 
resulted in questionable data since steady-state conditions : 
could not be maintained over the course of the individual 
data runs. 

During the early stages of this investigation, it was 
found that the steady-state conditions were fairly 
difficult to achieve. This was especially true after making 
a step change in the heat flux. For example, if the heat 
flux was increased, the cooling flow rate to the condenser 
must be increased to maintain the system pressure. Since 
the system showed about a 30-second time lag between the 
change of the valve position and the subsequent system 
response, the pressure was often seen to undergo cyclic 
variations. With considerable experience of the operator, 
it was possible to achieve steady-state conditions within 5 
to 10 minutes, as demonstrated by Reilly (Ref. 5). 

In order to Minimize the efforts necessary by the 
operator, two 600-W auxiliary cartridge heaters were 
installed within the boiling tee. With this arrangement, 


steady-state conditions could be maintained very easily by 
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maintaining the total heat load in the boiling tee ata 
constant value. For example, if the heat duty Pt o ps 


boiling tube was increased, the heat duty on the auxiliary 


heaters must be decreased by the same magnitude. іп this 
manner, steady-state conditions could be maintained with 
minimum attention. However, this arrangement required that 


the heat duty оп the auxiliary heaters be a maximum when 
the heat duty on the boiling tube is at a minimum. As 
determined by the early data runs (see Figure 2.6), the 
convective effects created by the bubbles of the auxiliary 
heaters created an artificially enhanced heat-transfer 


coefficient at low heat-flux settings. 


In order to alleviate the problem discussed 
above, while maintaining the advantages offered Бу 
the auxiliary heaters, the present investigator 
considered the use aL a convection shield. As 
shown in Figure DM this shield was positioned 
around the test tube, thus isolating it from the 


convective bubble effects created by the auxiliary heaters. 
The convection shield also provided increased 
mixing of the refrigerant and  condensate without 
interfering with the boiling phenomenon of the boiling 
tube or the circulation of the oil-refrigerant 
mixture in the boiling section by both the boiline 


tube and the auxiliary heaters. 
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р, DATA ACQUISITION AND REDUCTION 

The data-acquisition system used for this research was 
a Hewlett-Packard 9826A computer and a Hewlett-Packard 
3497A Automatic Data Асас 15111002” (012110 The НР-9826А 
computer was used to control the НР-3497А, to provide hard- 
copy storage capacity and to analyze the data obtained 
during this investigation. 

The data-acquisition system provided the means Бу 
which the researcher was able to review system performance 
and stability prior to the taking of any qar System 
parameters were monitored on the video screen following the 
scanning of all of the channels by the HP-S497A. Following 
the automatic scanning of all of the channels, the computer 
would display the desired saturation temperature (2.2 °C 
the observed saturation temperature, the liquid temperature 
at both ends of the boiling section and their average, the 
vapor temperature, and the temperature of the water- 
ethylene glycol mixture sump (-15 to -18 2С) Орсе 118 
steady-state conditions were met for the required satura- 
tion temperature within the prescribed variance (70.1 К) 
of the desired saturation temperature, the data-acquisition 
System Was initiated Гог the individual data run. 
Since there was a time lag in the request Тог the 
taking of data and the actual printing of the data due to 
the sampling of each channel twenty times before recording 


the data, the data-acquisition system would print the data, 
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allowing the researcher to review the results prior to 


storing the data. If, during the sampling of the data, a 


transient trend was observed due to "He cycling of the 
E12 refrigerant system or the temperature difference 
between the thermocouples measuring the liquid R-114 was 


not within the specified temperature range, the researcher 
had the capability to reject the affected data set. This 
capability allowed the researcher to apply а cursory 
analysis of each data point and ascertain i f the 
temperatures measured during that individual data set were 
within the prescribed limits of variance. 

The information scanned on each channel of the HP- 
3497A were thermal emf's provided by copper-constantan 


(type-T) thermocouples made of O.25-mm-diameter (30 gauge) 


wire. Two additional measurements were provided to the HP- 
3497A by the power-sensing device, as described by 
Karasabun  [Ref. 171, which converted the AC voltage dialed 


into the variac of the main boiling tube heater into a DC 
signal, allowing the HP-3497A to accurately scan the 
converted value and then allow the HP-9826A to 
automatically apply this value to the analysis calculation. 
Upon the initiation of the command to зге data, the 
data-acquisition system would sample each channel and then 
compute, according to the step-wise procedure outlined by 
Шагазабп (Ref. 17) and print the results. Appendix В 


contains a summary of the data-acquisition program anda 


37 


listing of the data-reduction program шиши! теа сит тон 


course of this research. Appendix C contains a 


representative sampie of the data printout received for 


each data run. 
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А. INSTALLATION OF TUBE IN BOILING APPARATUS 

The external surfaces of each of the boiling tubes 
were cleaned with Nitol (2 percent nitric acid and 98 
percent ethyl alcohol) solution and then йө with 
acetone to ensure tube cleanliness. This treatment was also 
effective in the removal of oil and grease residues from 
the tube surfaces. 

After air drying, the 101210 tube to be tested was 
installed in the boiling section of the apparatus. The 
system was sealed and then evacuated to approximately 29 
inHg by means of a portable mechanical vacuum pump (6), 
shown in Figure 2.1. The apparatus remained under the 
evacuated condition for a period of about one hour to check 
if there were any major leaks. System pressure was measured 


by means of a Marsh pressure gauge (30 inHg to 150 psi 


range with 5265 inHg and Et psi accuracy, 
respectively). If there was not a noticeable drop in the 
vacuum reading as observed on the pressure gauge, the 


system was charged with R-114 vapor, which was condensed Бу 
the  water-ethylene glycol mixture and then the absolute 
pressure was raised to approximately 0.186 МРа (27 psi). 


Upon reaching this pressure, an Automatic Halogen Leak 
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Detector, TIF 56067 with a sensitivity of 3 ppm 
concentration, was used to check for system leaks. 

Upon successful completion of the leak-detection tests 
and taking the necessary steps to fix any leaks found, the 
R-114 liquid level was adjusted to be 20 mm above the top 
of the boiling tube surface. Liquid-level adjustments were 


made by either boiling off the excess refrigerant to the 


pA reservoir (3), or drawing the required amount of 
refrigerant 11420394 antec the system from the R-114 
reservar The required . liquid level for boiling was 
maintained constant for all tubes tested by means of a 
permanent scribe mark on the inside of the oil intake 
endbell of the boiling section. At this initial refrigerant 
level, the system was now ready for boiling tests in pure 


refrigerant. 


B GENERAL SYSTEM OPERATION 

Table 2 provides a summary of all of the 83 data runs 
made during this investigation. Each data run is annotated 
to indicate the conditions involved in the respective run. 
The data runs were numbered sequentially with an indicator 
code for each tube preceding the number of the data run. 
Since this research is a follow up of previous research, 
the initial data run was indicated as number 168. All data 


runs were undertaken at a saturation temperature of 2/7250 
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КМ тепе о the boiling tube indicator codes for Table 2 


me as follows: 


WH -- Wieland Hard Copper Smooth Tube 
СТВ == Wieland Gewa-T (1.02 fins per mm) 
TXE -- Hitachi Thermoexcel-E 
THE -- Hitachi Thermoexcel-HE. 
Initially, all data runs were taken over a range of 


nine different heat-flux settings, ranging from 2 kW/m? to 
01 kW/m?. Additional heat fluxes of 0.5, 0.8, 1.2 kW/m? 
were later included for all tubes with enhanced surfaces to 
study the onset of nucleate boiling in the oil-refrigerant 
mixtures and to study the effects of hysteresis on these 
tubes. The initiation of all data runs was carried out in 
the same manner to ensure consistency in the obtained data. 
Each run was started by an initial system "warm-up" period 
allowing the R-12 refrigeration system to obtain a minimum 
coolant starting temperature of about -15°C. This was 
accomplished by lighting off of the R-12 system in 
conjunction with the effluent  recirculation/cooling pump 
(8) and slight opening of the contro! valve VC. This 


procedure allowed the system to gradually decrease system 


temperature and pressure, while minimizing the adverse 
effects of - shocking the boiling tube. The system was 
operated in this condition TTE a stable condition was 
реє: with the saturation temperature of 2.2°C. This 
Operation was most critical when the testing sequence 
involved increasing heat-flux conditions since 
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TABLE 2 
Summary of Data Runs 


Run % Heat No.of Data Auxiliary Convection 
Мо. Oi) Flux Points 2 Heater тэ 
WH168 О DEC 20) OTt without 
WH169 О INC 22 ce without 
WH170 1 DEC 20 off without 
WH171 1 I NC 22 off without 
WH172 2 DEC 29 OT withotm 
WH173 2 INC 20 ore without 
WH174 3 DEC 20 off without 
WH175 3 INC 22 orf without 
WH176 3 DEC 20 off without 
WH177 6 DEC 20 off without 
WH178 6 INC 21 of f withous 
WH179 10 DEC 18 of f without 
WH180 10 I NC MES off without 
WH181 10 I NC 12 Off у ії Є пен 
GTB182 О DEG 24 on with 
GTB183 О I NC 23 on with 
GTB184 О DEC 22 on with 
GTB185 1 DEC 29 on with 
GTB186 1 INC 22 on with 
GTBIST 2 DEC 34 on with 
GTB188 2 INC 18 on -with 
GTB189 © I NC 22 on with 
GTB190 3 DEC 20 on with 
GTB191 6 DEC 21 ог with 
GTB192 6 INC 20 on with 
GTB193 10 I NC 208 on with 
GTB194 10 DEC 21 on with 
WH195 O INC 22 on with 
WH196 О DEC 20 on with 
WH197 О DEC 17 оп with 
WH198 О INC 22 on with 
WH199 1 I NC 20 on with 
WH200 1 DEC 20 on with 
WH201 2 INC 20 on with 
WH202 2 DEC 20 on with 
WH203 3 INC . 29 on with 
WH204 3 TNC 20 on with 
WH205 3 DEC 20 on with 
WH206 3 I NC 21 on with 
WH207 6 DEC 20 on with 
WH208 6 I NC 20 on with 
WH209 10 INC 20 on with 
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TABLE 2 
Summary of Data Runs (cont'd) 


Run % Heat No.of Data Auxiliary Convection 
огт тэ Гол 5 Heater . | Shield 
WH210 10 DEC 21 on with 
TXE210 O DEC 18 on with 
ТХЕ211 О DEC 12 on with 
TXE212 O DEC T2 on with 
TXE213 O DEC 14 on with 
EXE214 О DEC 1272 оп with 
Пе 215 О I NC 20 on with 
TXE216 O DEC 14 on with 
ха не 1 INC 27 on with 
ШСЕ2І1З 1 DEC 14 on with 
TXE219 2 DEC 16 on with 
ТХЕ220 2 DEC 16 on with 
ТХЕ221 2 INC 25 on with 
TXE222 3 DEC 15 on with 
TXEZ23 3 INC 22 on with 
ера d I NC 2А on with 
ТХЕ225 6 I NC 30 on with 
ПКЕ 226 6 DEC 16 on with 
СЕТ 6 ІМС 25 on with 
DXD228 10 INC 24 on wort 
ТХЕ229 10 DEG 25 от with 
TXEZ230 10 INC 20 on with 
прса! O INC 20 on with 
THE232 O DEC 12 on with 
THE233 1 DEC 174 оп with 
ТНЕ234 1 INC 24 on with 
rt 235 1 ENG 27 on WIth 
THE236 2 DEC 14 on WIth 
EEIESS 2 INC 12 on with 
ТЕР 238 2 INC 8 on with 
THE240 3 DEC 14 on with 
THEZ241 3 INC 18 on with 
THE242 з INC 8 on with 
THE243 З ІМС S on with 
РНЕЗАА 3 INC = оп with 
пе 245 6 DEC 14 оп with 
THEZ246 6 I NC 24 on wa th 
ТНЕ247 6 I NC 18 on with 
ГПЕ2АВ TO DEC 14 on with 
THE249 10 I NC 24 on with 
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failure co Strict iy adhere to this voula result 135 
premature nucleate boiling and the subsequent purging of 
the data. This initial warm-up period required a minimum 
of one hour to ensure adequate system performance. 

Upon reaching stable conditions, the initial heat flux 
for the boiling tube was set by means of the resistor 
variaci Also, the heat flux for the auxiliary cartridge 
heaters was set by a second variac. The heat flux of the 
auxiliary cartridge heaters was limited to settings that 
allowed the heating of the liquid R-114 without creating 
excess turbulence on the boiling tube. As discussed 
earlier, the boiling tube was shielded from the effects of 


the auxiliary heater by means of the convection shield. As 


the Oa concentration of the refrigerant-oii mixture 
increased, the heat-flux setting for the auxiliary heater 
was additionally limited to lower heat fluxes. This 


precluded the generated oi! foam from carrying over the 
sides of the convection shield and interfering with the 
operation of the boiling tube. 


After setting of the desired heat flux for the data 


рал“ the system was monitored for the vapor temperature, 
sump temperature and the liquid refrigerant temperature 
(two readings; left side and right side of the boiling 
tube). Control valve VC was adjusted to allow the system to 
reach steady-state conditions. Additionally, the liquid 


refrigerant temperatures were compared to ensure minimum 
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variance. 
state conditions, 


to take data. 


The data-acguisition system 
channel, compute the 
the results for operator viewing. 
the prescribed limits, they were 
for permanent records. A 


printout is contained in Appendix B. 


After each run was 


the boiling section by means 


широ ГП valve VI. Table 3 


the required volume of ou 


proper mass percent of 2211 


puns. Mixing of the oil and 


accomplished by means of the 


ensure a homogeneous mixture. 


refrigerant mixtures was observed 


flux range of each data run, 


the foam over the convection 


highest heat flux. As the oil 


carryover increased, but was 


section of the apparatus. 


observed in the condenser 


Figures ed to 323 
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the data-acquisition system was 


would scan 


heat-transfer coefficient 


stored on 


completed, 
of 
contains 
added 


for each 


shield was 
concentration 
limited 

Carryover 
section 


demonstrate the 


ВЕ tilling the requirements of the steady- 


initiated 


each assigned 
7-8 (е ШИН 35:88 :8 


If the data were within 


a floppy disk 


representative sample of a data 


Oil was added to 


a graduated cylinder 


the listing of 


to obtain the 


series of data 


refrigerant liquid was 


auxiliary heater, to 


Foaming of the œil- 


over the entire heat- 


although carryover of 


limited to the 


increased, 


to the vertical 


oil was not 


of the apparatus. 


foaming action of 


the oll-refrigeramt mixtures at a heat flux ofi im 
kW/m? for oil concentrations of О, & and 10 Derc ent: 


Data runs were made for each boiling tube tested with 


pure R-114 and followed by runs of oil-refrigerant 
mixtures of increasing О content: Oil-refrigeroang 
mixtures öf 157208 3, 6 and 10 percent oil by mass were 
used for all tubes. 

TABLE 3 


OIL CONCENTRATIONS VALUES 


Percent Tota! Volume (cc) Step Change kce) 
Oil 00: Ve оте о мини 
Q -- «- = = 
1 26.9 2673 
2 54.2 2277 
3 82.2 28.0 
6 1659.6 87.4 
10 295. 2 129 16 
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Figure 53.1 


>N 
, 
5 
2 
x 
2 
£ 
£ 
5 

* 





Fo mine Action for Pure R-114 
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Enc 2 Foaming Action for R-114 with 
6 Percent Oil 
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IV. | RESULT 


—— —— 


AND DISCUSSION 


6-2 ықты 0 08 8088 0 —-—— "сю oe 


A. EFFECT OF CONVECTIONGSSEUNEISD 
Two sets of data runs were taken on the smooth copper 
tube; one set without the convection shield to validate the 


performance in comparison with the data obtained by Reilly 


(at a saturation temperature of -2.2°C) (Ref. 51; апа а 
second data run using the convection shield and the 
auxiliary cartridge heaters. These data runs allowed for 


an accurate comparison of the data taken over a period of 
time and provided justification for the use of the 
convection shield. As can be seen from Figures 4.1, ЖЕ 
ands, en 101005: З3 and 10 percent oil concentration, 
respectively, the present data with and without the shield 
in place are essentially the same. However, as discussed 
in section 11.С, the use of the shield provided consistent 
Steady-state conditions over the entire range of heat 
fluxes and oil concentrations, thereby minimizing the 
efforts required by the operator. 

Figure ze which contains the data тш for 
pure refrigerant for the smooth copper tube at a 
saturation temperature of 2.2°C, with and without the 
shield, demonstrates the repeatability of the data. 
Also plotted on this figure are Reilly's smooth-tube 


data for pure refrigerant at a saturation temperature of 
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D DE This data set represents a series of data 


runs, starting with а heat fru: Ron 100 kW/m and 
decreasing the heat flux to 3a minimum of 1 kW/m?. The 
three sets of curves are nearly identical with the 


deviation between the Reilly data and the most recent data 
being attributed to the pressure effect and the possible 
aging of the tinned interface between the thermocouple 
sleeve and the inner surface of the smooth copper tube. 
The maximum discrepancy of the computed wall superheat 
between Reilly's data and the present data is about twenty 
percent at the maximum heat-flux setting. This figure also 
Shows the typical uncertainty of about 10 percent computed 


for the wall superheat values during this investigation. 


B BOILING PERFORMANCE TOF THE SHOOT ВЕ 
Figure 4.4 shows the performance of the smooth copper 


tube in pure R-114 and refrigerant-oil mixtures of 3 and 


10 mass percent oil. The behavior of the smooth copper 
tube in pure К-114 liquid demonstrates typical 
nucleate pool-boiling performance. As can be seen in 
Figure 4.4, the line between points A and B represents 
a region of constant slope (оп а log-log scale) for 
the variation of heat flux with the tube wall superheat 
(Two-Tsat). This section of the curve represents natural 
convection  (i.e., no bubbles were generated). 
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Upon reaching point B (see Figure 4.4). the boiling 


tube wall superheat begins to decrease even though the 
meat flux is continually increased. At this point, the 
tube begins КО demonstrate the characteristics of 
incipient nucleate boiling. This region is known as the 
mixed boiling region, where the transition from natural 
@emvection to nucleate boiling occurs. This event is 
characterized by the ` increasing number of active 
nucleation sites with increasing heat flux. It should be 
noted that during all data runs, the unheated ends of the 
tubes did not show any nucleation sites. Reilly CRef. 5] 
also Stated that the unenhanced ends of the boiling tube, 
at au observed heat fluxes, underwent оштру natural 
convection despite a smal! amount of heat axially conducted 
along the tube walls. It was observed that this transition 


luum natural convection to nucleate boiling occurs very 
quickly, usually within a few seconds after the activation 
of the first nucleation site. 


At the heat flux represented by point C in Figure 


4.4, all the available activation sites of the boiling tube 
appear to be active and the wall superheat begins to 
increase with increasing heat flux. From рит ~ to 


point D, where the heat flux is increased to a maximum 
ої 91 kW/m-, new activation sites are generated (maximum 


heat flux was limited by cartridge heater rated output). 
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As the oil concentration was increased, there was a 
significant increase in the amount of foam generated. 
122 obtaining a maximum heat flux at point D, the 
heat flux was decreased and is С as the region 
between points D and E on Figure 4.4. This region of the 
curve represents decreasing wall superheat; however, more 
activation sites that were created at the higher heat-flux 
settings remain active to lower heat fluxes. This 
increased range over which the boiling tube continues to 
nucleate results in an increase in the heat-transfer 
performance since the wall superheat temperature for the 
nucleate boiling at any particular heat flux is less than 
the tube wall superheat for the natural convection heating 
at the "S heat flux. 


~ 


Stephan [КеЁ. 9] reported that the effects of adding 


eu introduces a mass-diffusion resistance to the liquid 
and therefore lowers the heat-transfer coefficient. Reilly 
(Ref. 5l, in his discussion of his smooth-tube data at a 


Saturation temperature of -2 212780 1 ё гс 4.5) “stated that 
the 3 mass percent oil and 10 mass percent oil 
concentrations demonstrated higher wall superheat 


temperature than that observed in the pure 8-114 liquid. 


However, Reilly observed that the wall superheat for the 3 
percent oil concentration was higher than that for the 
10 percent oil concentration. He attributed this to the 
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non-linear physical property characteristics of the 


metrigerant-oil mixtures. As can be seen in Figure 4.4, 
the general trend of the data is the same with the 
exception of the 3 percent Gil -concemeration wal | 
superheat, which nearly overlays the data for pure 8К-114 
mi guid. This deviation from the Reilly data may be due in 
part to the difference in the saturation temperatures. The 
sddaition of oil to the pure R-114 liquid, up to a maximum 


Of 10 percent, should delay the transition from natura! 
convection to nucleate pool boiling on the tube. With a 
mefrigerant-oil mixture of 10 mass percent oil, it was 
observed that the surface nucleation sites did appear ta 
spread more slowly as the heat flux was increased. This 
condition continued until the tube was fully engulfed in 
the nucleate boiling mode. This observation was also Е Т 
By Reilly Ке 5] апа иа ггееѕ міт the contention of 
Chongrungreong and Sauer (Ref. 8] and Thome (Ref. 6] that, 
with the exception of the surface tension, the non-linear 
variations in the physical properties of the refrigerant- 
oil mixtures explain the variation in the boiling heat 
meamsfer of the tube. 


Figure 4.6 shows the heat-transfer coefficient of the 


smooth tube in refrigerant-oil mixtures, plotted as a 
Шс топ of heat flux. Figure 4.7 shows the same plot for 
Reilly's ааба Тот а saturation temperature of  -2.2?C 


59 


OT 


'equ] u3ooug лој з3чәтотјјәо) лој5џелј-зеон Зитттоц 


Cau/M» b 


„О? 





HLOOWS 


HLOOWS 


HLOOWS 


зак | 
вал]. 


9'y eunSrg 


OT 
OT 


xn 
О 
= 
OI 
E = 
№ 
6 
Ж 
ы 
ова 01 Z 
oag Є O 
oer О ~ 
Xn t4 110 
зоэн "© 
OT 


60 


OT 


ES 


Сә ЧИМ? 


b 


TO 
pro 
Ir 
150 
110 
IFO 


“2әч1 езеа әчпі чцзоошс 5,Л41112еч 
JO] BPUSFOTJJaoD Јәјѕиејлј-зеән Вит ТОН 


quadiad 
yuəo1iəd 
диззлод 
yuəoiəd 
quadiad 


зи20184 


°з еаабта 


owe Co 


ОТ 
9 


t 
с 
І 
0 


c] uu. зро ыы -.- 








is) 


CH p/m? 


БЯ 


[Rete Sri As can be seen, the two graphs are nearly 
identical, representing 2556 repeatability of the data. 
Also, the effect of the addition orf up to © percent Tons 
(Figure 4252) consistently represents an approximate 


decrease of 10 percent tn the heat-transfer coefficient. 


As the oil concentration was increased from 6 со 10 
percent, the decrease in the heat-transfer coefficient was 
a maximum of about 35 percent. These decreases іп the 


heat-transfer coefficients were dependent on the comparison 


of the pure К-114 and the 10 percent oil concentration 
at equal heat-flux settings. 
22 BOILING PERFORMANCE OF THE СЕМА-Т TUBE 


Figure 4.8 shows the nucleate pool-boiling performance 
of the GEWA-T (1.02 fins/mm) tube in R-114-oi1l mixtures 
The magnitude of the tube-vwall superheat was considerably 
less than the values obtained for the smooth tube discussed 
itm. (Section ув This lower wall superheat сап be 
attributed to the specially formed fins of the GEWA-T tube. 
The performance of the GEWA-T tube is similar to that of 
the smooth tube with the curves of the increasing heat flux 
to the point of incipient nucleate boiling (point В 
paralleling the performance of the smooth copper tube 
(Figure 4.4). The same condition exists throughout the 


entire range of decreasing heat fluxes. 
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The region bounded by points A and B shows the 
characteristic  natural-convection heating behavior. The 


onset of incipient nucleate boiling at point B is clearly 


defined for the pure К-114 liquid. However, as the oil 
concentration was increased, the point of incipient 
nucleate boiling becomes less clearly defined and the 


region of mixed boiling (point B to С) increased with 
[im@ereased heat flux. It was also observed that in the 
region from A to B, the 3 percent oil mixture showed the 
smallest wall superheat and the pure refrigerant showed the 
highest wall superheat. The reason for this behavior is 
mee fully known at present, but may be due to convection 
effects within the channels of the GEWA-T tube which are 
enhanced by the foaming action in the presence of oil. 

The region нан by the points C and D for the 
increasing heat-flux condition is indicative of fully 
developed nucleate pool boiling. The region between points 
D and E represents nucleate pool boiling for decreasing 
heat Trux: Figure 4.8 shows that, as the oil 
concentration was increased, the wall superheat increased 
at a given heat flux. This is contrary to the results for 
the condition of increasing heat flux in the region of 
natural-convection heating. 

When comparing the performance of the smooth tube 


"igure 4.4) to that of the GEWA-Tube (Вог аа). the 
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maximum adverse effect of the increased oil content for the 
smooth tube, occurs at the higher heat-flux settings. For 
the GEWA-T tube, the maximum decrease in the tube 
performance for increased oil content, occurs in the lower 
-and intermediate heat-flux range. 

Figure 4.9 shows the heat-transfer coefficient as a 
function of heat flux over the range of  R-114-oil 
mixtures (Сал m 2 3, 6 and 10 percent oil). Data obse Sd 
@ver the range of R-114-o01! mixtures of 1, 2, 3 and 6 
percent oil are nearly coincident. The 10 percent oil 
concentration data represent an approximate 30 percent 
reduction of the heat-transfer coefficient when compared to 
the heat-transfer coefficient obtained during the data for 
pure R-114. Figure 4.10 shows the direct comparison of the 


GEWA-T tube with the smooth tube for R-114-oil mixtures of 


©, dj and 10 percent oil. The GEWA-T tube represents an 
overall enhancement of 4.2 over the baseline smooth copper 
tube at a heat flux of 40 kW/m? with pure R-114 and an 
overall enhancement of 3.8 with 1O percent oil. These 
results compare favorably with the results of Yilmaz, 


Hwalek and Westwater (Ref. 18) for a 12.3-mm OD GEWA-T tube 
in p-xylene, in which they obtained an enhancement of 5.3 
in the pool-boiling heat-transfer coefficient over their 
baseline smooth copper tube. Additional ly, Yilmaz and 


Westwater [Ref. 12] tested a similar surface in isopropyl 
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alcohol in which an enhancement of 2.0 was observed. Wantage 
and Lepere (Ref. 13] tested а 17.9-mm OD GEWA-T surface in 
both R-113 and FC-72 and reported enhancementsmof 2359000 
2.95, respectively. 

Figure 4.11 shows more clearly the degradation 
that oil causes in the boiling heat-transtem 
performance for the СЕМА-Т tube. This figure pir 


the heat-transfer coefficient ОТ the GEWA-T tube in 


the ca mixture relative to the heat-transfer 
coefficient 1n риге R-114 as a fune tian of oil 
concentration. From the shape of the curves, the 
oil-caused degradation does not show a clear trend 
with the oil concentration at all practical heat 
fluxes. Over the range of oil concentrations, except 
at a heat flux setting of 37 КЫ/т?, there is only 
a small degradation (up to a maximum of 35%) in the 
boiling performance. This is probably due to the 
large channels and reentrant cavities in the tube, 
allowing the "pumping" action of the bubbles to contu" 
tO remove the 'oril-rrch liquid throughout the ем 
heat-flux range without а ‘Significame reduction mim 
performance. The only Signi f reane decrease in 
performance occurred at heat fluxes greater than 5 


kW/m? with. 10° percento 
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D. BOILING PERFORMANCE OF THE FTHERDMOEXCED.POBUSID 


Figure 4.12 shows the nucleate роо І -Бот Үе 
performance of the Thermoexcel-E tube in R-114-0 
mixtures. As can be seen in the figure, the region of 


natural-convection heating for the condition of incressmm 
heat flux, points A to B. for the pure R-114 is very та 
with the point of incipient nucleate boiling Occurring ean 
reasonably low heat flux (1.5 kW/m?). АХ this роги ШОШ 
incipient) nucleate Dol ling; the tube performs in the mixed 
boiling region as a limited number of nucleation sites 
become -active. From ‘point Corto Ci; the tube appears to 
Operate in more of a natural-convection heating mode as no 
new nucleation sites are activated, although previously 
activated sites continue to nucleate. This phenomenon is 
probably due to the inability of the cartridge heaters to 
produce uniform heat flux at these lower settings. At 
point G with a heat-flux setting ОТ 5 kW/m:, more 
nucleation sites are generated and as the heat flux is 
increased to 60 kW/m?, the tube becomes fully nucleated. 

Por. the 3 <“pereene,oin concentration the point Мєн 
incipient  nucleate boiling, point B', decreases as the non- 
linear physical properties of the R-114-oil mixture allow 
the onset of the mixed boiling region to occur. Also, the 


region of the mixed boiling, point B: to D P is 
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far larger than that for pure R-114,; as the generatione 


new nucleation sites is retarded by the increased oil 


content: 

For the 10 percent oil concentration, on the обе 
hands the point of incipient nucleate boiling, point B GSE 
is delayed considerably. At approximately 14 kW/m?, the 
tramsition from natural-convective heating to nucleate 


boiling occurs. The exact reason for the difference between 
the 3 percent and 10 percent behavior is not known. however 
it is possibly due to the change in the physical properties 
of the R-114-072] mixture. 


The region of nucleate pool boiling for the condition 


of decreasing heat flux represented by points D їо Е. р! 
тоң FES and D'' toT E for О, З and 10 percent oil 
concentrations, respectively, for the Thermoexcel-E tube, 


represent classic nucleate pool-boiling characteristics. 


AS the oil concentration increases, the wall superheat also 
increases throughout the range of heat-flux settings. At 
the lower heat fluxes, the largest increase in the wall 
superheat occurred at 3 percent oil concentration when 


compared to the pure R-114, while the 10 percent oil 
concentration was observed to undergo continual degradation 
in tube performance over the entire range of  heat-flux 
settings, although a more noticeable reduction occurred at 


the higher heat fluxes. At a heat-flux setting of 91 


T2 


kW/m?, the R-114-oil mixture of 10 percent oil was observed 
to undergo a reduction in tube performance by a factor of 
з, For the practical heat-flux range of ЗО to 40 kW/m?, 
there is an increase in wall superheat by factors of 1.4 
ис Бо for. R-114-oil mixtures of 3 and@ 10 percent, 
respectively. 

Figure 4.13 shows the heat-transfer coefficient of the 
Thermoexcel-E tube іп refrigerant-oil mixtures, plotted as 
a function of heat flux. As can be seen in the figure, the 
heat-transfer coefficient remains nearly constant 
throughout the range of heat-flux settings for pure R-114, 
while there are slight reductions for the 1, 2 аан: 
percent oil concentrations. Additionally, it can be seen 
that at these lower oil concentrations, the heat-transfer 
coefficients are nearly coincident. As the 8-114-01! 
mixture oil concentration is increased to 6 percent, the 


tube performance increases with increasing heat flux until 


reaching a heat-flux setting of 60 kW/m?, at which time 
tube performance begims to deteriorate. This is contrary 
ШШ Пе tube performance for O, 2, and З percent oil 


concentrations in which the heat-transfer coefficient was 
observed to increase slightly over the range of heat 
fluxes. Form che 10 percent -oil concentration, the heat- 
transfer coefficient decreased as the heat flux was 


increased. 
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Figure 4.14 shows а direct comparison of the 
Thermoexcel-E tube to the smooth copper tube for R-114-oil 
"E tures of O, 3 and 10 percent. The Thermoexcel-E tube 
presents enhancements of 6.6 and 4.1 ver the performance 
of the smooth copper tube in pure R-114 and an R-114-oil 
mixture of 10 percent oil, respectively, at a heat-flux 
setting of 40 kW/m?. The results for the pure R-114 compare 
favorably with the results obtained by Yilmaz, Palen and 
Taborek (Ref. 19] for the Thermoexce!-E in p-xylene. 

Figure 4.15 shows the degradation that oil causes in 
the boiling heat-transfer performance for the Thermoexcel-E 
tube. This figure plots the ratio of the heat-transfer 
coefficient of the Thermoexcel-E tube with oil to the value 


without oil as a function of oil concentration. The shapes 


of the curves continue to demonstrate the non-linear 
variation of the physical properties of refrigerant-oi! 
mixtures. However, the significant degradation of the tube 
performance is probably due to the inability of the 


generated bubbles to scavenge the interior tunnels of the 
Thermoexcel-E tube. Ihis IS іп contrast to the СЕМА-Т 
tube, where more scavenging is possible due to its larger 


channel size. 


Е. BOILING PERFORMANCE OF THE THERMOEXCEL-HE TUBE 
Figure 4.16 shows the nucleate BOG lS bor] time 
performance of the Thermoexce!-HE tube in R=-114-0il 
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mixtures. In the region between рроїпЕз Ага сг 006 
natural-convection heating occurs for increasing heat ГІШ 


Although similar in performance to the Thermoexcel-E tube, 


the ТРегмоехсе | -НЕ tube presented a more linear 
eharacteérisStric than that of the Thermoexcel-E tube, 
representing а more uniform heat distribution. At point m 
on. she guess the point of incipient nucleate boiling for 
pure R-114 liquid occurs at approximately 2 kW/m?. FOT 
R-114-07 | mixture of З апда 10 percent oil concentrations 
спе ТӘСІ oF incipient nucleate boiling, point B'. осећа 
at approximately 8 kW/m?. Also in this region, the pure 


К-114 and the Ае mixture gro з percent oil 
concentration are nearly coincident in the range of heat 
flux from 500 W7m* to 2 kW2me. 


From роте В tor c; the Thermoexcel-HE tube in the 


pure R-114 liquid achieved full nucleate boiling at 
approximately 37 kW/m?, with intermediate data points in 
the region of mixed boiling. These intermediate data 


points indicate that the Thermoexcel-HE tube was not able 


to activate all nucleation sites simultaneously for 
increasing heat-flux conditions. Тһе ЮК-114-о11 піхічге 6% 
10 percent oil was able to achieve activation of all 


nucleation sites upon reaching 22 kW/m?, without any data 
points contained in the mixed boiling region. With XS 


percent oil, the Thermoexcel-HE tube was able to achieve 
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ТТМ nucleate botling at 37-7 mee The Тћегтоехсе не 


tube 111 the R-114-oil mixtures appeared to activate 


all nucleation sites simultaneously. 

The region from points C to D and f r om 6 ОР 
represent the nucleate pool boiling for increasing heat 
f Tux. The performance of the Thermoexcel-HE tube in this 
region appears to decrease as the oil content was 
increased. The 10 percent oil concentration increase цн 
wall superheat by a factor of approximately 2 above the 
wal | superheat for the pure R-114 liquid for heat-flux 


settings of 40 kW/m?. 

The regions between the points D and E and D' апа: 
represent nucleate pool-boiling for decreasing heat-flux 
conditions. As can be seen from the figure, the greatest 
reduction in performance occurs in the lower heat-f lux 


settings between the pure R-114 and the 3 percent oil 


concentrstron. This is similar to the results гат the 
Thermoexcel-E tube, where the tube performance decreased 
as oil content increased. There is a significant decrease 


in the Thermoexcel-HE tube performance over the range of - 
heat-flux settings as the oil concentration is increased to 
10 percent. 

Figure 4.17 shows the heat-transfer coefficient of “ШЕ 
Thermoexcel -HE tube in refrigerant-oil mixtures, plotted 
as a function of heat mt I ux. As can be seen, the 


performance 15 very similar to that of the Thermoexcel-E 
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збе» (Figure 4.13). As was observed for the Thermoexcel-E 
tube, the Thermoexcel-HE tube undergoes a slight reduction 


in the heat-transfer coefficient as the oil concentration 


was increased from 1 to 3 percent. As the oil 
concentration is increased to 6 апа 10 percent oil, the 
performance reduction is quite noticeable. With the 10 
percent oi! concentration, tube performance decreases at 
both the lower and higher heat-flux settings. This, s 


contrary to the bindings for the Thermoexce!-E tube where 
the performance decreased over the entire range of heat- 
Dux settings. from low heat flux to high heat-flux 
settings. 

Figure 4.18 shows a comparison of the Thermoexcel-HE 
tube to the smooth copper tube at R-114-oil mixtures of O, 
З апа 10 percent oil concentrations. The Thermoexcel-HE 
tube displayed an enhancement of 6.1, 4.6 and 4.1 at a heat 
flux setting of 40 kW/m? and R-114-oil mixtures of O. 3 and 
10 percent oil, respectively, over the same parameters for 
the smooth copper tube. These results are nearly identical 
to those found for the Thermoexcel-E tube. 

Figure 4.19 shows the degradation that the oil has on 
the boiling heat-transfer performance of the Thermoexcel-HE 
tube. At the lower heat fluxes, the effect of the oil 
concentration is slightly more pronounced than that 


observed for the Thermoexcel-E tube throughout the range of 
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Oil concentrations, although more severe at the lower oil 
concentrations. At the higher heat fluxes, the trend 


appears to reverse. 


СО COMPARISON TO REILLY'S DATA FOR THE HIGH FLUX TUBE 


Figures 4.20, 4.21 and 4.22 provide a comparison of 


the smooth, High Ёл Хо, GEWA-T. Thermoexcel-E and 
Thermoexcel-HE tubes for O EE апа 10 percent oil 
concentrations, respectively. For pure R-114 liquid, the 


High Flux,Thermoexcel-E and Thermoexcel-HE tubes operating 


at a heat flux of approximately 20 kW/m?, show nearly 
identical performance. As the heat flux was increased, the 
performance of the High Ех tube continued to 
increase while the performance of  Thermoexcel-E and 
Thermoexcel-HE tubes remained nearly constant. At heat- 


flux settings lower than 20 kW/m?, the Тһегпоехсе!-Е апа 
Thermoexcel-HE tubes showed higher heat-transfer 
coefficients than that for the High Flux tube. As the oil 
concentration was increased to 3 percent, the performance 
of all five tubes decreased, as can be seen in Figure 4.21. 
ал іопаї 1у, the GEWA-T tube's performance appeared to 


decrease the least Since, as discussed earlier in section 


КЛ С, the finned surface of the GEWA-T tube is not 
adversely effected by the small increase in the oil content 
due to its large reentrant cavities. From this 
observation, it can be stated that the performance of the 
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tubes with porous coatings or smal! reentrant cavities is 
adversely affected by the increase in oil content due to 


the generation of an oil-rich film in the nucleation sites 


impeding the tube's performance. This is further evidenced 
pCigure 4.22. which depicts the performance of all five 
tubes їп ЗЛА = оти mixtures of 10 percent oil. The 
performance of all the tubes 15 reduced significantly, 


while the performance of the High Flux tube is severely 
decreased in the range of heat-flux settings from 22 kW/m? 
to the maximum setting of 91 kW/m?. At approximately 45 
kW/m?, the performance of the High Flux tube, the GEWA-T, 
the Thermoexcel-E and the  Тпегтоехсе! НЕ tubes are 
approximately equal, although the performance of the 
Thermoexcel-E and the 3 Thermoexcel-HE tubes remain nearly 
constant over the range of heat-flux settings of 14 kW/m? 
to 50 kW/m*. Although the performance of the High Flux 
tube exceeds that of the GEWA-T, the Thermoexcel-E and the 
Thermoexcel-HE tubes over this range, its performance 
begins to severely degrade at the midpoint of the heat-f lux 
range indicated. In the heat-flux range above 45 kW/m?, 
the GEWA-T tube outperformed all the other tubes. This is 
attributed to the ability of bubbles generated in its 
larger channels to continually scavenge the oil-rich layer 
from the tube, and therefore remove the insulating oil from 


the tube's channels and exterior surface. 
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A. 


V. CONCLUSIONS AND RECOMMENDATIONS 


————— - v —— зе па — и он смиље - — — n 


CONCLUSIONS 


In pure R-114, at a heat flux of 30 kW/m? the 


pool-boiling heat-transfer coefficient of the 
Thermoexcel-E, Thermoexcel-HE and GEWA-T tubes 
were approximately ое, 6.1 and 4.2 times, 
respectively, larger than the baseline smooth tube. 
Ina R-114-0il mixtures of 3 percent oil, at a 
heat flux of ЗО kW/m?, the pool-boiling heal 
transfer coefficient of the Thermoexcel-E, Thermo- 
excel-HE and GEWA-T tubes were approximately 4.8, 
4.6 and 4.0 times, respectively, larger than the 


smooth tube under the same conditions. 


In R-2114-01] mixtures of 10 percent oil, at a heat 
flux of 30 kW/m?, the performance of the Thermoexcel- 
E and Thermoexcel -HE tubes were nearly identical, 
with an enhancement of 4.1 over the pool-boiling 
heat-transfer coefficient of the smooth tube. The 
GEWA-T tube, under the same conditions, produced an 
enhancement of 3.8 over the smooth tube. 


As oil content in the refrigerant-oil mixture 
increased, the performance of enhanced surfaces 
decreased; however, the amount of this decrease was 
directly related to the size of the subsurface 
channels (1.e., the larger the reentrant cavity ог 
channel, the less adverse effect of the increased oil 


content). 


As the oil content increases in the channels, the 
greater the adverse effect of the oil-rich layer has 
on the surface, insulating the surface from the 


refrigerant and therefore decreasing the performance 
of the surface. 


The GEWA-T surface in refrigerant-oil mixtures, with 
its larger subsurface channels, was able to 
scavenge the oil-rich layer from its cavities by 


the bubbling action generated within the channels. 


The performance of the Thermoexcel-E and Thermoexcel- 
HE tubes, with their smaller pore size and smaller 
tunnels, showed significant performance decreases 
as oil content was increased. The smaller pore and 
tunnel size under these conditions, were unable to 
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effectively disperse the oil-rich film from the sub 
-surface cavities and therefore allowing the oil-rich 
film to insulate the tube's surface. 


RECOMMENDATIONS 

The single-tube testing should be expanded to include 
multi-tube bundles. This would provide more realistic 
data with regard to the effects of tubes operating in 


close proximity to each other. 


A uniform heat flux is necessary within the boiling 


tube and within the boiling section itself. There- 
fore, the present cartridge heaters need to be 
replaced with more precise heaters, capable of pro- 


ducing a uniform heat flux in the axial direction. 


The testing of surfaces should be expanded to study 
the effects of various refrigerant oils. 


Conduct experiments to establish the physical proper- 
ties of refrigerant-oil mixtures to enable further 
understanding into the boiling mechanisms involved. 


There needs to be additional studies undertaken on 


the effects of additional boiling temperatures on the 
mixture properties of the refrigerant-oil mixtures. 
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APPENDIX A 


THERMOCOUPLE CALIBRAT [ON tr eto 


Karasabun [Rete 171 described the thermocouple 
calibration procedúre in вова се = mm Reilly (Ref. 2 


stated that the data-reduction program utilized differences 


in thermocouple readings im all computations, such as the 
amount oof the tube wall superhest (Two: = паи Еш Reillv 
stated "that tne calibration of the thermocouples would be 


necessary if the thermocouples were measuring quantities 


such as absolute temperature to determine actual fluid 
properties, The procedure utilized by Karasabun and Reilly 
TOT the initial instrumentatoone or the experimental 


apparatus is as contained in Appendix A of Reference 5. 

To determine the calibration for the thermocoup vec! 
two thermocouples were made, one from the beginning of the 
wire roll and the second from the end. The seventh-oraer 
polynomial equation provided by the wire manufacturer was 
corrected by adding to it a second-order curve fit of the 
Variation of the manufacturer's computed temperature for a 
given emf from a known set of reference temperatures as 
measured by a Hewlett-Packard 2804A quartz thermometer. The 
thermometer had а temperature resolution of =0.0001 K and 


ат ассыгасу оғ 8а” 
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The manufacturer's  emf to temperature conversion 
equation is: 
fms Sone tf As es tae Roe AM EMT UN (А.1) 


asES + а, Е? + а; Е’ 


where 

E = thermocouple reading in volts 
T = temperature 
ao = 0.100860910 
a, = 25727.94369 
a2 = -767345.8295 
pou" 75025595 361 
aa = -9247486589 
as =  6.97688x10'! 
РГА = 
8: - 2.94078х10:% 


Figure A.1 (CRef. 5] shows the difference in the quartz 


thermocouple readings and the actual thermocoupie readinzs 


(discrepancy) plotted as a function of temperature. AS can 
be seen from the figure, the two thermocouples agreed to 
within 205208 of each other. The manufacturer's 


calibration equation needed an increase of approximately 


21 К to more accurately convert the emf's to the true 
temperature. The second-order polynomial correction was as 
follows; 

ВЕБ DS t DT + bT? (AL 2) 
where 
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DCP = discrepancy (K) 


bo =  85.62589568« 107 
b. =. 3. 7619902107 
Б: - -5 059925205053 
T? = thermocouple readings from 
equation Awl In C 
Therefore, the temperatures computed by the data-reduction 


program СрЕР5, were emf's converted directi to 
temperatures by equation A.1 and corrected by equation A.2. 
THis calibration was. thus estimated to result m 
temperature readings within 20405 K ай the true 


temperature. 
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The 


APPEND. B 


DATA=REDUCTION PROGRAM [REF. 5) 


data-reduction program below consists of the 


following sections: 


Main Program - Menu of subprogram options 

Sub Main - Takes or reprocesses data 

Sub plot - Plots data on log-log scale 

РОТУ - Computes least squares curve fit 
Su D or ein - Plots data on linear-linear scale 
Sub Stats - Computes average and standard 


deviation of data 


Subprogram Main consists of the following steps: 


17 


2. 


Creates data file for data and plotting file. 
Annotates tube type. 

Monitors heat flux or saturation temperature to 
establish steady-state conditions. 

Scans HP3497 channels and stores in data file. 
Converts raw emf's to temperatures, current and 
voltages. 

Computes the heat-transfer rate for the tube's 


internal cartridge heater. 


Computes the average wall temperature of the 
tube, the average wall superheat and 
the film temperature. 


96 


10. 


Dis 


eZ 


3. 


Computes the physical properties of the 
R-114 using the given correlation temperature. 
Computes the natural-convection heat-transfer 
coefficient of R-114 for the non-boiling 
tube ends. 

Computes the heat loss from the non-boiling 
tube ends. 

Calculates the corrected heat fluxes from 
the testing tube to the R-114. 

Calculates the boiling heat-transfer coefficient 
of the R-114 from the tube being tested. 

Prints the data and stores the heat flux and tube 


wall superheat values in the plot file. 


The following is the listing of the complete data-reduction 


program 


(DRP5) written in Basic 3.0 for the Hewlett-Packard 


9826 computer. 


S 


1000! 
1003! 
1006! 
1009! 
1012 
1015 
1018 
1021 
1024 
1027 
1030 
1033 
1036 
7055 
1042 
1045 
1048 
1051 
1054 


EtCIS) 


1057 
1060 
1063 
1066! 
1069 
1072 
1075 
1078 
1081 
1084 
1087 
1030 
1093 


— СОА) © 9) 00 O м Б — OO UM 


Бо О) ММГ) ~. — 


FILE NAME: DRPS 
DATE: Qerober 19221334 
REVISED: April 25. 1986 


COM /Idp/ Idp 

PRINTER IS 1 

CALL Select 

INPUT “WANT TO SELECT ANOTHER OPTION (1=Y,0=N)?",Isel 
eo THEN GOTO 1018 

B 


BEER 
PRINTER Slo. | 
Nae “DATA COLLECTION/REPROCESSING COMPLETED" 


SUB Main 

СОМ /1Іар/ Іар 

COM ЕСА ЕСЕ ТЕАТ 

COR 7Wil7y D2.Di,Do.L Lu Ken 
DIM Emf(C12),TC12) ,Dfat112).D2a(112).Diía( 11), Doa(C 1177. саст) Еа) Кеса а 
,Ї п5(4)(15) 

DATA 0.10086091,25727.94368.-767345.8295.78025595.81 
DATA -9247486589,6.97688bE*11,-2.66192E*13,3.94078E *14 
READ С(») 

DATA "Smooth"."High Flux"*,"Thermoexel-E","Thermoexel-HE" 
DATA Smooth,Hiah Flux, Turbo-B,Hiagh Flux Mod.Turbo-B Mod 
READ Tnš(*) 

PRINTER IS 701 

BEEP 


IF Idp=4 THEN PRINTER IS 1 

IF Idp=4 THEN GOTO 1774 

INPUT "ENTER MONTH. DATE AND TIME ш DD:HH:MM:SS)",Dates 
OUTPUT 709:°TD*:Dates 

OUTPUT 709;"TD* 

ENTER 709: Dates 


PRINT 
PRINT " Month. date and time :"';Dates$ 
PRINT 
ES USING "10X,"*MOTE: Program name : DRP4*"*"* 


INPUT "ENTER DISK NUMBER” ,Dn 

EDI USING "16X.""Disk numper = "",22":0п 

BEEP 

ey "ENTER INPUT MODE (0=#3054A.1=FILE)*,Im 
INPUT "SELECT HEATING MODE (O=ELEC: 1!=WATER)*,Ihm 


INPUT “ENTER THERMOCOQUPLE TYPE (O=NEW,1=O0LD)*,Ical 
IF Im=0 THEN 

BEEP 

INPUT "GIVE A NAME FOR THE RAM DATA FILE".D2 files 
PRINT USING *16X,""Neu file пате: "",148":02_# | 1е5 
512е1=20 

CREATE BDAT D2 files$.Sizel 

ASSIGN 9Р:1е2 TO 02_+ 1125 


DUMMY FILE UNTIL Nrun KNOWN 
01 files2"DUMMY" 

CREATE BDAT D! files.Sizel 
ASSIGN «File! TO D1 files 
OUTPUT ЗР:1е1:Оасез 
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1270 
1273 
1276 
1213 
1282 
1285 
1288 
1291 
1294 
1237 


IF Ihm*0 THEN 

BEEP 

INPUT "ENTER NUMBER OF DEFECTIVE TCS (OsDEFAULT)D",Idtc 
IF Idtc*0 THEN 

Ldtc!=0 

Ldtc2=0 

PRINT USING "16X,""No defective TCs exist'"" 

END IF | 

IF Idtc-1 THEN 

BEEP 

INPUT *ENTER DEFECTIVE TC LOCATION" ,Ldtc! 

PRINT USING "16X.""TC is defective at location "",D'";Ldtc!l 
Ldtc2=0 

END IF 

IF Idtc=2 THEN 

BEEP 

INPUT “ENTER DEFECTIVE TC LOCATIONS*.Ldtcl!.Ldtc2 

PRINT USING "16X,""TC are detective at locations "",D.,4X,D"sLdtct,Ldtc2 
END IF 

IF Idtc>2 THEN 

ВЕЕР 

PRINTER IS 1 

BEEP 

PRINT “INVALID ENTRY” 

PRINTER IS 701 

GOTO 1180 

END IF 

END IF 

OUTPUT @Filel;Ldtcl1,Ldtc2 


| Im*! option 


ELSE 

BEEP 

INPUT “GIVE THE NAME OF THE EXISTING DATA FILE",D2 fileS 
PRINT USING "16Х,“”014 file name: °",14A":D2_files 

ASSIGN *File2 TO D2 files 

ENTER 8Fi1le2:Nrun 

ENTER 3F11e2:DoldS$ 

PRINT USING "16X.""This data set taken on : "",148";DoldS 
ЕШПЕН SETle2:ULdtcl.Ldtc2 

IF Ldtc!>0 OR Ldtc2>0 THEN 

PRINT USING "16X.""Thermocouples uere detective at locations:"",2(3D,4X)0"; 


elite | ЕЈ 


1300 
1303 
1306 
1309! 
1312 
1315 
1318 
1321 
1324 
1327 
1330 
1333 
1336 
1359 
1342 
1345 
1348 
1351 


END IF 

ENTER &SFile2;Itt 

END IF 

IF Im-0 AND Ihm-1! THEN 15935 
BEEP 


INPUT “WANT TO CREATE A PLOT FILE? (02N,17#Y)*,Iplot 
IF Iplot=1 THEN 

BEEP 

INPUT “GIVE NAME FOR PLOT FILE",P files$ 
CREATE BDAT P d uH 4 

ASSIGN 8Plot TO P files 

END IF 

IF Ihms! THEN 

BEEP 

INPUT “WANT TO CREATE Uo FILE? (0-N, 1*0 ",Iuf 
23 Bute 1 THEN 


EE 
INPUT “ENTER Uo FILE NAME" .Ufiles 
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1354 
1357 
1360 
135: 
1366 
1559 
(272 
1375 
1378 
1381 
1384 
1387 
1390 
1993 
1396 
1399 
1402 
1405 
1408 
1411 
1414 
1417 
1420 
1423 
1426 
1429 
1432 
1435 
1438 
1441 
1444 
1447 
1450 
1453 
1456 
1459 
1462 
1465 
1468 
1471 
1474 
1477 
1480 
1483 
1486 
1489 
1492 
1495! 
1498! 
1501! 
1504! 
1507 
1510 
1513 
1516 
1519! 
1522! 
1525 
1528 
1531 


CREATE BDAT Utiles.4 
ASSIGN 
END IF 


BEEP 
INPUT 


ВЕЕР 
INPUT 


END 
END 


PRINT 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


INPUT 


IE 

Ir 
PRINTER IS 1 
IF Im-Q THEN 
BEEP 


"file TO Ufiles$ 


"WANT TO CREATE Re FILE? (0=N.1=¥)",Ire 
IF Ire=! THEN 


"ENTER Re FILE NAME" ,Refiles 
CREATE BDAT Refiled.10 
ASSIGN ®Refile TO Refiles 


USING 


USING ' 


USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 
USING 


USING 
USING 
USING 
US ING 
USING 


Itt 


*áX, “Select tube number" `° 
IF Ihm=0 THEN 


OUTPUT $Filel;Itt 


END IF 
PRINTER 
IF Itt«10 THEN PRINT USING *18X,""Tubpe Number: 
IF Itt»9 THEN PRINT USING "16Х, * "Tube Number: 


IF Ihm-1 THEN PRINT USING "16X,"*Tube Type: 


ВЕЕР 
INPUT 
ВЕЕР 


ELECTRIC HEATED MODE 


IS 701 


Smooth 4 inch Ref"*"* 

Smooth 4 inch Cu (Press/Slide)*"* 
Soft Solder 4 nen Goes 

Soft Solder 4 inch HIGH FLUX°°” 
Wieland Hard 8 inch"*" 

HIGH FLUX 8 inch""* 

GEHA-X 19 Fins/in"*"** 

GEHA-K 25 Fins/in"'""* 

GEHA-T 19 Fins/in"** 

GEHA-T 26 Fins/in"'** 
THERMOEXCEL -E"** 
THERMOEXCEL -HE" ** 


Smooth tube""" 
High Flux""" 
Тигро-В""" 
High Flux Mod*"*"* 
Тигбо-8 Мод“"" 


ИВЕ 
~ 002510: 
79, 547 2 ne Cie 


"ENTER QUTPUT VERSION (0s-LONG. 1! =5Н0ЕТ ,2=М0МЕ) “ , Том 
INPUT “SELECT (0=LIQ.1*VAP.2=(LIQ+VAP)/2)".Ilav 
DIMENSIONS FO TESTED TUBES 


Di=Diameter at thermocouple Боз 1005 


DATA .0111125. 


.Ü111125..0111125) 


0129540..012446..0129540, .0100965 


рата .0100965,.01157, 01157) ‚01157. . 01157 
READ Dia(*) 
D1=Dia( Itt) 


D2=Diameter of test section to the base of 


DATA .015875, .0153875. 015375, 


fin 
.015824. 015375 "015824. .01270 


DATA 0127, -0138,.0138..0138..0128 
READ р2а(») 
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1534! 

1537! Di*Inside diameter of unennanced ends 

гэ тег. Ољам..0127,.0132..0127,.0132. .0111125,.0111125,..0118,.0113..01 
18..0118 

1543 READ Dia(*) 

1546! 

1549! DosQutside diameter of unenhanced ends 

НО 11 015925. .015875..015875,.015824,.015875..015824,.01270,.01270,.01331, 
нозі .01331..01331 

1555 READ Doa(*) 

1558! 

1561! L*Lenath of enhanced surface 

Нева о бАГО .1016..1016,.1016,.1016,.2032. 2032. 2032. 20532. 2002, 200 Нара ка 


32 

1567 READ La(*) 

1970! 

1573! Lu-Lenath of unenhanced surface at the ends 

IBS 35878 ..0254,.0254,.0254,.0254,.0762..0752,.0762,.0762,.0762..0762,.0762,.07 


62 

1579 READ Lua(*) 

1582! 

1585! Kcu-lhermal Conductivity of tube 

1588 DATA 398,344,344,45,344,45,344,944,398,398,398, 398 
1591 READ Со 

1594 [Е Ihm=1 THEN 

15537 ! 

1600! Data statements for water heating mode 

1603! 

2:5 7958/8 0.019875,0.015875.0.0169.0.0138.0.0169,0.0.0,0.0 
1609 READ D2a(*) 

По DAIA 0.0127.0.0127,0.0145,0.0127,0.0145,0,0.0.0.0 
1615 READ Dial») | 

I рата 0.015875.0.015875,0.0169.0.015875.0.0159,0.0.0,0,0 
1621 READ Doa(*) 

1624 DATA 0.3048,0.3048,0.3048,0.3048,0.3048.0,0,0,0.0 
1527 READ (а(») 

16530 ОҺТА 0.0254,0.0254,0.0254.0.0254,0.0254,0,0,0.0.0 
1623 READ Lua(*) 

СЕБ DATA 398.45,398,45,398.0.0,0,0.0 

1639 READ Kcua(*) 

1642 END IF 

1645 02=02а(1+%) 

1648 Di-Dia(Itt) 

1651 DosDoa(Itt) 

(654 Le*La(Itt) 

1657 Lu-Lua(Itt) 

1660 Kcu-Kcua(Itt) 

1663 Ха”.8 

1666 Fr*.3 

1669 IF Itt-0 THEN Cf=!.70E+9 

1672 IF Itt>0 THEN Cf-3.7037E*10 

1675 Rñn=PI*(Do `2-Di 2)/4 

1673  P-PI*Do 

1681 IF _Ihm=1 THEN 

1684 BEEP 

1687 INPUT “TUBE INITIATION MODE. (1=HOT HATER.2=STERM,3=COLD HATER)",Itim 
1630 IF Itims! THEN PRINT USING "16X.""Tube Initiate: Hot Hater""" 

1693 IF Itim-2 THEM PRINT USING "16X.""Tube Initiate: Steam'"" 

1696 If Itim=3 THEN PRINT USING "16X.""Tube Initiate: Cold MWater""" 

1699 INPUT “TEMP/VEL MODE: (0*T-CONST.V-DEC:!*T-DEC,V-CONST; 2=T-INC,Y-CONST)"., 
Itv 
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IF Itvsó THEM PRINT USING "15X.""Tsmp/Vel Mode: Constant/Decreasinqg"'"" 
IF Itv=1 THEN PRINT USING “16X,.""Temp/Vel Mode: Decreasing/Constant"”" 
IF Itv22 THEN PRINT USING "168X.""Temp/Vel Mode: Increasina/Constant""" 
INPUT “WANT TO RUN MILSON PLOT? (12Y,0-ND",Ivil 

IF Ihm=1 AND Iwil=0 THEN 

IF Itt-0 THEN Ci-.032 

IF Itt- OR Itte3 THEN Cr=.059 

IF Itt=2 OR [tt=4 THEN Ci=.062 

Beer 

INPUT “ENTER CI (DEF: WH=.032,HF=.059,TB=.062)",Ci 

PRINT USING "TBX. "Sieder-Iste 777 

PRINT USING "16X,"" Constant з "”,2.40 081 


END IF 


END IF 

IF Ihm-! AND Im-2! AND Iuii1-! THEN 
IF Itt=0 THEN Ci=.032 

IF Itt-1 OR Itt-3 THEN Ci-.059 

IF Itt-2 OR Itt-4 THEN Ci-.062 
ASSIGN SFile2 TO + 

CALL Wilson(Cf,Ci) 

ASSIGN @File2 TO D2 file$ 

ENTER 3Fi1le2;Nrun,DoldS.Ldtc!,Ldtc2,Itt 
END IF 

Nsub=0 

IF Idp=4 THEN 1һт=! 

IF Ihm=1 THEN Nsubd=8 


Nte=6 

IF Ihm-0 THEN Ntc-!2 

341 

Sx=0 

Sy =0 

Sxs=0 

Sxy=0 Я : 


Repeat: ! 


IF In*0 THEN 

Dtld-2.22 

Іаоч2 

ON KEY 0.15 RECOVER 1801 

PRINTER IS 1 

PRINT USING "4X.""SELECT QPTION""* 

PRINT USING "S5X.""0-TAKE DAaTa""" 

IF Ihm=0 THEN PRINT USING *"6X.""1-SET HEAT FLUX""" 
IF Ihm=1 THEN PRINT USING *6X."*"1=SET WATER FLOW НАТЕ""" 
PRINT USING@GA."“2=Ser Їєаэ 277 

PRINT USING "4X.""NOTE: KEY 0 - ESCAPE""" 

ВЕЕР 

[МРОТ Ido 

IF Ido»2 THEN Ido-2 

IF Ido-0 THEN 2296 


! LOOP TO SET HEAT FLUX OR FLOWMETER SETTING 


IF Ido-! THEN 

IF Ihm=0 THEN 

OUTPUT 709:“АВ АЕ12 85:13 VRS” 

ВЕЕР 

INPUT “ENTER DESIRED Odo” ,Dadp 

PRINT USING "4X,""DESIRED Qdp ACTUAL Qdp""" 
Err-1!000 

FOR I-1 TO 2 

QUTPUT 709:"85 Sa" 
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1882 
1885 
1888 
1831 
1834 
1897 
1900 
1903 
1906 
1309 
1912 
1915 
1518 
1921 
1924 
[927 
1930 
1933 
1936 
1939 
1942 
1945 
1948 
1351 
1954 
1957 
1960 
1963 
1966 
1969 
[872 
1975 
1978 
1981 
1984 
1887! 
1990! 
|593 
1455 
1599 
2002 
2005! 
2008 
2011 
2014 
2017 
2020 
2023 
2026 
2123 


бит= 0) 

ҒОВ Јі=1 Т0 5 

ENTER 709:E 

Sum*Sum-*E 

NEXT Ji 

ЇЕ 121 THEN Volt=Sum*5 
IF I=2 THEN Amp=E 

NEXT I 

Aadp=Volt+#Amp/ (PI*D2*L) 
ТЕ ABS(Aqdp-Dadp)>Err THEN 
IF Aqde>Dadp THEN 

BEEP 4000,.2 

BEEP 4000,.2 

BEEP 4000,.2 


EN 

PRINT USING "4X.MZ.3DE,2X,MZ. 3DE" :Dadp, Aadp 
WAIT 2 

GOTO 1876 

ELSE 


BEEP 

PRINT USING °4X,MZ.3DE,2X,MZ.3DE”";Dadp,Aaqdp 
Err=500 

WAIT 2 

GOTO 1876 

END IF 

ELSE 


ВЕЕР 

INPUT "ENTER FLOWMETER SETTING",Fms 
GOTO 1819 3 
END IF 

END IF 


LOOP TO SET Tsat 

IF Ido#2 THEN 

IF Ikdt={ THEN 2011 

BEEP 

INPUT “ENTER DESIRED Tsat",Dtld 

PRINT USING “4X,"" DTsat ATsat Rate Ty Rate""" 
Ikdt=1 

014140 

014240 

Мп=1 

Nrs=Nn MOD 15 

Мп=М№п +1 

IF Nrs={ THEN 

PRINT USING "4X.** Tsat Trat Tld2 Tv Tsump 


тонг" ~ 


2032 
2035 


2044 
2047 
2050 
2053 
2056 


END IF 

IF Ihm=9 THEN OUTPUT 709:“AR AFS AL1I1 YRS" 
IF Ihm=1 THEN OUTPUT 709:“AR AFO ALS VRS" 
РОК 1=1 ТО 5 

IF Ihm=0 AND I>4 THEN 2125 

биті) 

QUTPUT 709:5А5 58" 

ЕОК Ла1 10 20 

ENTER 709:Elia 
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Tinlet 


Tpile 


2170 
21125 


2113 
2182 
2185 
2188 
21:! 
2134 
2194 
2200 
2203 
2206 
2209 
22 
2215 
2219 
2221 
2224 
2227* 


Sum=SumrE lia 

NEXT Ji 

Eliq=Sum/20 
TLd=FNTvsv(Eliq) 

IF I=1 THEN Tldl2Tld 
ЇЕ [=2 THEN Tld2-2Tld 
IF [=3 THEN Tv=Tld 

IF [#4 THEN Tsump-Tld 
IF I*S THEN Tinlet=Tld 
IF [#6 THEN Tout=Tld 
NEXT I 

IF Ihm=t THEN 

OUTPUT 709:“AR AF20 AL20 VRS” 
OUTPUT 709:"AS SA” 
Sum=() 

FOR Kk=1 TO 20 

ENTER 709;E 

Sum=SumtE 

NEXT Kk 

Emt (7) =ABS(Sum/20) 
Terle=Emt (7)/3.96E-4 
END IF 
Atld=(Tld1+Fld2)«.5 

IF ABS(AtlLd-Dtld)>.2 THEN 
IF Atild>Dtld THEN 

ВЕЕР 4000, .2 

ВЕЕР 4000, .2 

ВЕЕР 4000,.2 


ВЕЕР 250..2 
BEEP 250. .2 
ВЕЕР 250..2 
END IF 
Еггізб6%14-0141 
01414-4614 
Ёгг242-Ту-0142 
Old2=Tv 


IF [hm=0 THEN PRINT USING *4X.S(MDD.DD.2X)"*:Dtld.Tld1.Tld2.Iv.Tsumo 

IF Ihm=1 AND Idp=0 THEN PRINT USING °4X,7(MDD.0D,2X)“3:Dtld.Tld1,JTld2.JTv.Ts 
ump.Tiniet.Tpile 
2176* IF Ihm-! AND Idp-4 THEN PRINT USING "4X,5(MDD.DD.,2:0,3(M30.DD, 2X0 " 5Dtld. TI 
dl,Tld2,Tv.Tsump,. Tin 


WAIT 2 
GOTO 2020 
ELSE 
IF ABS(Atld-Dtld)>.1 THEN 
IF Atld>Dtld THEN 
BEEP 300052 
BEEP 3000..2 
БЕЗЕ 
BEEP 300..2 
BEEP 800,.2 
ШЕ 


Errt=Atld-Old!1 
0191 =а4+1а 
Err2*Tv-0l1d2 
Old2=Tv 


IF [hm=0 THEN PRINT USING "4X.S(HDD.DD.2X)";Dtld.T1d1,T1d2.Tv.,Tsump 
IF Ihm=1 THEN PRINT USING "4X,5(MDD.DD,2:09,3(M3D.DD. 1:0 "^ ;Dtl1d. Tldl.Tl1a2,. Tv 
.Isump,Tinlet,Tpile, 
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э нг em, eae 


2230 
223 

2236 
2239 
2242 
2245 
2248 
2251 
2254 
2057 


WAIT 2 
GOTO 2020 
ESE 


BEEP 
Ёгг14-0:14-0141 
Oldt=Atld 


0192=Т\ 
IF Ihm-0 THEN PRINT USING "4X,5(MDD.DD,2X0 ":;Dtld,Tldt,Tl1a2.Tv, Tsumo 
IF Ihm=1 THEN PRINT USING *4X,8(MDD.DD,2X)";Dtld,Tld1,Tld2, Tv, Tsump,Tinlet 


»Tpile, Tout 


2260 
2263 
2266 
2269 
2272 
2275! 
2278 
2281 
2284 
2287 
2230! 
2233! 
2296 
2239 
2302 
2505 
2308 
2311 
2314 
2317 
2320 
232: 
2325 
2423 
2:32 
2:135 
2138 
2341 
2344 
2347 
2350 
2553 
2356 
2859 
2362 
2-55 
2368 
2371 
2374 
СС 
2380 
2383 
2386 
2389 
2399 
2335 
2338 
2401 
2404 
2407 


WAIT 2 

0010 2020 

END IF 

END IF 

END IF 

ERROR TRAP FOR Ido QUT OF BOUNDS 
IF Ido>2 THEN 

BEEP 

GOTO 1819 

END IF 


TAKE DATA IF Im=0 LOOP 
ГЕ Ikol*1 THEN 2308 


BEEP 
шад ENTER BULK OIL Вор 
ol= 
IF Iha-0 THEN OUTPUT 709:°AR AFO ALI1 VRS" 
IF Ihm=1 THEN OUTPUT 709;3"AR AF0 ALS VRS” 
IF Ihm=0 THEN Ntc=!2 
FOR I=1 TO Nte 
QUTPUT 709:°AS SA” 
Sum=0 | 
РОВ Jis! TQ 20 Е 
ENTER 709:E 
Sum*Sum*E 
IF I=(9-Nsub) OR [2(10-Nsub)° THEN Et(Ji-!)*E 
NEXT Ji 
Kd1=0 
IF I-(9-Nsub) OR Is(10-Nsub) THEN 
Eave=Sum/20 
Зит=д. 
FOR Jk=0 TQ 19 
IF ABS(Et(Jk)-Eave)<S.0E-6 THEN 
Sum=SumtE t( Jk ) 
ELSE 


Ка! =Ка1 +1 

END IF 

NEXT Jk 

IF I-(9-Nsub) OR I=(10-Nsub) THEN PRINT USING "4X,""Kdl » "",DD";Kdl 
IF Kdl>10 THEN 

BEEP 

BEEP 

PRINT USING "4X.""Too much scattering in data - repeat data set""" 
GOTO 1816 

END IF 

END IF 

Emt (I) =Sum/(20-Kd1) 

NEXT I 

IF Ihm*! THEN 

QUTPUT 709;"AR АҒ20 8120 VRS" 
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2531 
2393 
2536 
2559 
2502 
2505 
2608 
2611 
2614 
2617 
2620 
2623 
2626 
2623 
2632 
263 
2638 
2641 
2644 
2647 
2650 
2653 
2656 
2653 
2662 
2665! 
2668! 
26/1 
2674 
26/7 
2680 
2683 
2686 
2683 
2632 
2595 
2638 
2701 
2704 
Y 
2710! 
204131 
2115 
2113 
2722 
2725 
2/28 
2791 
2/34 
2737! 
2/40! 


NEXT I 
Tu*Tya/(83-Idtc) 
END IF 
Tld-T(9-Nsub) 
Tld2=T(10-Nsub) 
Tldaz(Tld+Tld2)+*.5 
Ту«Т(11-МецЬь) 

IF Itt«3 АМО Ihm-0 THEN 
Tlg29-99.99 

Ім= (7010) +7011) )/2 
END IF 
Tsump*T(12-Nsub) 
IF Ihm=0 THEN 2635 
Tinlet*T(13-Nsub) 
Tout*T(14-Nsub) 

IF Ihm=0 THEN 
Amp=ABS( Ir) 

Чо 1 с-АВ5 (4 г) »25 
Q=Volt*Amp 

END IF 

IF Itt=) AND Ihm=0 THEN 
Kcu*FNKcu(Tu) 

ELSE 

Kcu*Kcua(Itt) 

END IF 


FOURIER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
IF Ihm=0 THEN Tus-Tu-Q*L0GCD2/D12/(2*PI*Kcu*L) 
IF Ilav-0 THEN Tsat=Tlda 

IF Ilqv=! THEN Tsat=(TldatTv)*.5 

IF Ilav22 THEN Tsat=Tv 

IF Ihm=! THEN 

Тама-Тіпіеє . 
Grad=37.9853+.104388*Tava 

Тагор=ав  (Ет! (7))»1.2+6/(10»бгад) 
Tavacs*Tinlet-Tdrop*.5 

IF aBS(Tava-Tavqc)»5.01 THEN 
Тама-(ТаматТамас)».5 

GOTO 2689 

END IF : 


COMPUTE WATER PROPERTIES 
IF Ihm=1 THEN 

Kw=F NKw( Tava) 
Миша=Е ММиш (Тама) 
Срш=РМСрш( Тама) 

Ртш= ҒМРтш( Тама) 

КРош=Е МВРош ( Тама) 
Тш: = Тама 


Compute MDOT 


2 ,dot-3.9637E-3*Fms* (S. 6 T95SE-3-r ms (B. B2DU BESTE 
219) 


27461 
2749 
2752 
2155 
2758 
2761 
2764 


g ot dot*(1.0965-Tinletz(1,35pa E In S O S P 
t= 

Q-Mdot*Cpu*Tdroo 
LmtdeTdrop/LOG((Tinlet-Tsat)/(Tinlet-Tdrop-Tsat)) 
Uo*Q/(PI*0o*L *Untd) 

Ru*Üo*LOG(Do/Di?/(2. *Kcu) 

Tus[sat*Fr*U mtd 
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2410 
2412 
2416 
2419 
2422 
2425 
2428 
2431 
2434 
2437 
2440 
2443 


2578 
2581 
2534 
2587 


OUTPUT AgS:"AS Sa" 

Sum=0 

FOR Kk=! TO 20 

ENTER 709;£ 

Sum*Sum*E 

NEXT Kk 

Emt (7) =ABS (Sum) /20 

END IF 

IF Ihm-0 THEN 

OUTPUT 709: "АК АҒ12 Аі 13 VRS" 

FOR I=! TO 2 

OUTPUT 709;“А5 SA" 

Sum=0 

FOR Ji=1 TO 2 

ENTER 709:E 

Sum=Sum+€ 

NEXT Ji 

IF I=1 THEN Vr=Sum/2 

IF I*2 THEN Ir=Sum/2 

NEXT I 

END IF 

ELSE 

IF Ihm=0 THEN ENTER 8F:ile2:Bop,ToldS,Emt(*) Чг, Іг 

T or THEN ENTER @File2;3;Bop.ToldS,Emt(*) ,Fms 
DI 


CONVERT emt’S TO TEMP,VOLT ,CURRENT 
Twa 


SFOR l=! TQ Ntc 


IF Idtc»0 THEN 

IF I-Ldtc! OR I-sLdtc2 THEN 
T(I)=-99.99 

GOTO 2536 _ 
END IF 

END IF 

IF Itt<4 AND Ihm=0 THEN 

IF I>4 AND I«9 THEN 

T(I) =-99.99 

GOTO 2536 

END IF 

END IF 

TCI) =FNTvsv(Emt (I) ) 

NEXT I 

IF Itt«4 THEN 

FOR I=1 TO 4 

IF I-Ldtc! OR I-sLdtc2 THEN 
Тша=Тша 

ELSE 

Тша"Тша%Т(Т) 

END IF 

NEXT I 

Тша=Тша/ (4-1Іаїс) 

ELSE 

IF Ihm-! THEN 2599 

FOR I=1 TO 8 

IF I=Ldtci OR I=Ldtc2 THEN 
Тша = Тша 

ELSE 

Тша= Тша+Т (1) 

END IF 
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2257 
2220 
27213 


272729 
2782 


2805 
2809 


2815 
2819 
2821 
2824 
2827 
2830 
2933 
2836 
2833! 
2842! 
2845 
2848 
2851 
2854 
2857 
2860 


2866 
2869 


VusMdot/(Rhou*PI*01'2/4) 
Rew=Rhow*Vw*01/Muwa 
HisCi*Ku/Di*Reu ' .3*Pru'^(1/2.) *(Muua/FNMuu( Tu1))^.14 
Tuic-Tavgq-Q/(PI*Do*L *H1) 
IF ABS(Twi-Twic)>.01 THEN 
Tuis(Tui*Tuic)*.S 

GOTO 2773 

END IF 

Tuis(Tui*Tuic)*.5 
Hos1/(1/Uo-Do/(Di*Hi)-Ru) 
END IF 

END IF 

IF Ihm=1 THEN 
Thetab2Q/(Ho*xPI*Do*L) 
Tw*Tsat+Thetab 


Thetab=Tu-Tsat 
END IF 

IF Thetab«0 THEN 
BEEP 


INPUT *"THALL«TSAT (O0-CONTINUE. !-END)",Iev 
IF lev=0 THEN GOTO 1804 

IF Тем»! ІНЕМ 3130 

END IF 


COMPUTE VARIOUS PROPERTIES 
Tfilm=(TuwtTsat)*.5 
Rho#FNRho(Tfilm) 
Mu=FNMu( TF ilm) 
K=FNK(Tfilm) 
Co*FNCo(Tf ilm) 
Beta=FNBeta(Tfilm) 
Hfg=FNHfg( Tsat) 
Ni=Mu/Rho 

Аірһазк/ (КһожСр) 
Pr*Ni/Rlpha 
Psat=FNPsat(Tsat) 


COMPUTE NATURAL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNENHANCED END(S) 

Hbar*!90 

Fes(Hbar*P/(Kcu*0Q))^.S5*Lu 

Tanh*FNTanh(Fe) 

Thetas*Thetabp*Tanh/Fe 

Xx2(9.81*B8eta*Thetabp*Ü0o 3*Tann/(Fe«Ni*01pha))^.166667 
Ту=(1+(,559/Р7) (9/16)) (8/27) 
Hbarc*K/Do*(.6*.387*Xx/Yy) '2 

IF aBS((Hbar-Hbarc)/Hbarc)».001 THEN 
Hbar=(Hbar+Hbarc)*.S 

GOTO 2890 

END IF 


! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 


Ql-(Hbar*P*Kcu*0)^ .S5*Thetab*Tanh 
Qc*0-2*01 
Aas-PI*D2*L 


| COMPUTE ACTUAL HEAT FLUX AND BOILING COEFFICIENT 


QOdp=0c/As 
Htube=Qdp/Thetab 
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248) Csfs(Cp*Thetab/Hfq)/(Qdp/(Mu*Htq)*(.014/(9.81*Rho)/.5)^(1/3.) *Pr 1,7) 
2950! 

2953! RECORD TIME OF DATA TAKING 

2956 IF Im=0 THEN 

2959 OUTPUT 709;"TD* 

2962 ENTER 709: ToldS 

2965 END IF 

2968! 

2371! OUTPUT DATA TO PRINTER 

2974 PRINTER IS 701 

2977 IF Iov-0 THEN 

2980 PRINT 

2983 PRINT USING "10X.""Data Set Number * "",DDD.2X.""Bulk Oil X = "",р0.0,5Х,1 
48" : Ј,Вор . Т0145 

2986 IF Ihm-0 THEN 

2989 PRINT USING "10X,""*TC No: 1 2 S 4 5 6 7 


EE SRI USING "10X,"*Temp :5",8(1Х,МрО.003:1(1),Т(2),1(3),Т(4),Т(5),Т(6),Т6 
) 


2995 PRINT USING "10X,"" Tua таа Tiied2 Tsapr Psat Tsump""* 

Б РКІМТ USING *10X,2(MDD.DD, 130, 1X, MDD. DD. 1X,2( TX, MDD.DD) ,2X,MDD.D* s Tu, T1d,T 

ld2,Tv.Psat.Tsump 

3001 PRINT USING "10X,"" Thetab Htube Шар” 

3004 PRINT USING "10X,MDD.3D,1X.MZ.3DE, 1X. MZ.3DE" 3 Thetab,Htube,Qdp 

3007 ELSE 

3010 Bu USING "10X,"" Fms Vu Tsat Tinl Tdroe Thetab q Uo 
mos" 


3013* PRINT USING *'10X,4(2D.DD.130 .Z.3D, !1X,.DD.DD. 1X, 3CMZ.3DE, 10 * s F ms, Vu, Tsat, Tai 
nlet,Tdrop.Thetab,Qd 

3016 END IF 

3019 END IF 

3022 IF Iov-1 THEN 

3025 IF J=1 THEN 5 

3028 PRINT 

3031 IF Ihm=0 THEN 

3034 PRINT USING "10Х, “* КОМ No 0:1% Tsat Htube Оде Тһеғавһ""” 

4037 ELSE 

3040 PRINT USING "10X."" FMS OILZ TSAT HTUBE QDP ТНЕТАВ" ®" 
3043 END IF 

3046 END IF 

3049 IF Ihm=0 THEN 

222 ERNUS ING =12X,3D0,.4X,DD.2X.MDD.DD.3(1X,MZ.3DE)":J.Bop,Tsat,Htube,Qdp,The 
5 

3055 ELSE 

m PRINT USING "12X.3D.,4X.DD,2X,MDD.DD,3C1X, MZ.3DE? ";Fms, Boe, Tsat,Htube,Qdp.T 
eta 

3061 END IF 

3064 END IF 

3067 IF Im-0 THEN 

3070  BEEP 

3073 INPUT "OK TO STORE THIS рата SET (12Y,0-2N)2?*.,0k 

3076 END IF 

3079 IF Ok=! OR Im=! THEN J=J+} 

3082 IF Ok=1 AND Im=#0 THEN 

3085 IF Ihm=0 THEN OUTPUT ЗЕ: Іеї: Вор, То! 4$, Ет? (+), Чг, Іг 

3088 IF Ihm-! THEN OUTPUT ӨЕ: іеї ; Вор, То! 4$, Еп? (+) ,Етѕ 

3091 END IF 

3094 IF Iuf=! THEN OUTPUT 3Ufi1l1e:Vu.Uo 

3097 IF [re=1 THEN OUTPUT GSRefile;fms.Reu 

3100 IF (Im=! OR Ok=!) AND Iplot=! THEN QUTPUT @PLlot:Qdp,Thetab 

3103 IF Im=0 THEN 


2199 


plot 
3199 
3202 
3205 
3208 
3211 
3214 
3217 
3220 
3223 
3226! 
3229! 
3232 
3235! 
3238 
3241 
3244 
3247 
3250 
3252! 
3256 
3253 
3262 
3265 
3268 
32711 
3274 
3277 


ВЕЕР 

INPUT “WILL THERE BE ANOTHER RUN (127,02N)?",Go_on 
Nrun2J 

IF Go опе-д THEN 3130 

IF Go on«»0 THEN Repeat 

EESE 

IF J«Nrun*! THEM Repeat 

END IF 

IF Im*0 THEN 

BEEP 
ү USING *"10X,""NOTE: "".Z2Z."* data runs were stored in file "DOM = 
iles 

ASSIGN sFile! TO * 

OUTPUT ФЕ: е2: Мгип- 1! 

ASSIGN ЭЕ:1е! ТО 01 бі|еб 

ENTER 83Fi1lel:Date$S.Ldtcl.Ldtc2.Itt 
OUTPUT F ile2;DateS,Ldtcl,Ldtc2,Itt 
FOR I=1 TO Nrun-1 

IF Ihm=0 THEN 

ENTER 3File!i;S8op,ToldS.Emt (*),Vr,Ir 
OUTPUT 3F11e2;Bop,Told$,Emt C*) Ут, [г 
ELSE 

ENTER 3File!:8op.ToldS.Emf (€) ,Fms 
OUTPUT $Fi11e2;Bop,ToldS,Emf e) ,Fmns 
END IF 

NEXT I 

ASSIGN File! TO = 

PURGE “DUMMY” 

END IF 

BEEP 

PRINT 

IF Iplot=! THEN PRINT USING "10X.""NOTE: "*,Z2Z."" X-Y paırs were stored ın 
data file "",1 : 
ASSIGN @File2 TO + 

ASSIGN *Plot TO =» 

IF Iuf=! THEN ASSIGN S3Ufile TO =» 

IF Ire=1 THEN ASSIGN ®Refile TO + 
CALL Stats 

EEP 

INPUT “LIKE TO PLOT DATA (1=Y,0=N)?7",O0k 
IF Око! THEN CALL Plot 

SUBEND 

CURVE FITS OF PROPERTY FUNCTIONS 

DEF FNKcu(T) 

OFHC COPPER 250 TO 300 K 

Tk=T+273.15 IC TO K 

K*434-.112*Tk 

RETURN K 

FNEND 

DEF FNMu(T) 

170 TO 360 K CURVE FIT OF VISCQUSITY 
Tk =T+273.15 IC TO K 

Ми=ЕХР (-4.4636+(1011.47/ТК) ) *1.0Е-3 
RETURN Mu 

FNEND 

DEF FNCo(T) 

180 TO 400 K CURVE FIT 909 Ср 

ТК =1+273.15 ІС ТО К 

Ср=.40188+1. mo J3*Tk*1.51494E-6*Tk ^2-6.67853E- 1O* Tk ^3 
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3280 Ср-Ср»1000 

3283 RETURN Co 

3296  FNEND 

3289 DEF FNRho(T) 

3292 Tk-2T*273.15 IC OTO K 
3295 X=1-(1.8*Tk/7S3.9S) !K TO R 
3238 Ro=368.32+51.148414*X*(1/3)+16.41801S*X+17.478838*X ° °.S+1.119828*X 2 
3301 Ro=Ro/.062428 

3304 RETURN Ro 

3307  FNEND 

3310 DEF FNPri(T) 

3313 Pr=FNCp(T)*FNMu( T)/FNK(T) 
3316 RETURN Pr 

3319 FNEND 

3322 DEF FNK(T) 

3325! T«360 K WITH T INC 

3328 K=.071-.000261+T 

3331 RETURN K Е 

3334  FNEND 

3337 DEF FNTanh(X) 

3340 P=EXP(X) 

3343 Q=1/P 

3346  Tanhz(P-9)/(P*Q) 

3349 RETURN Tanh | 

3352  FNEND 

3355 DEF FNTvsv(V) 

эг COM /Ce/ С(7),Їса! 

3361 T=C(0) 

3364 FOR I-! TO 7 

3367 T=T+C(I)*V I 

3370 NEXT I 

3373 IF Ical=1 THEM 

3376 Т-Т-6.7422934Е-2%Т%(9.0277043Е-3-Тж%(-3.3259917Е-5)) 


ELSE 
3392 TsI*8.526897E-2*T*(3.76199E-3-T45.0689259E -5) 
3385 END IF 
3388 RETURN T 
3391  FNEND 
3394 DEF FNBeta(T) 
3397 Rop=FNRho(T+.1) 
3400 Rom=FNRho(T-. 1) 
3403 Beta=-2/(Rop+tRom) *(Rop-Rom)/.2 
3406 RETURN Beta 
3409 FNEND 
3412 ПЕҒ ҒМНҒа(Т) 
3415 ]Hfq*!.3741344E*5-T*(3.2094361E*2*T*1.2165143) 
3418 RETURN Hfg 
3421 FNEND 
. 3424 DEF FNPsatiTc) 
3427! 0 TO 80 deg F CURVE FIT OF Psat 
3430 Tf-1.8*Tc*432 
3433 Pa-25.945525*Tf*(,15352082*Tf£*(1.4840963E-3*T£*9.6150671E-95)) 
3436 Ра-Ра-14.7 
3439 IF Pq>0 THEN | *2PSIG.-2in Ha 
3442  PsatzPq 
3445 ELSE 
3448 Psat=Pq*29.92/14.7 
3451 END IF 
3454 RETURN Psat 
3457  FNEND 
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3460 
3463 
3466 


DEF FNHsmooth(X.Bop.Isat) 

DIM ACS) ,.B(S) ,C(S) .D(5) 

DATA .205265..25322..319048. .55322. .73909.1.00258 
рата .74515..729892. .73189,. 71225, .68472. .64197 
рата .41092..17726..25142.,.54806..819816,1.0845 
рата .71403, .72313. . 72555. .6365631, .555867 ,. 51889 
READ AC*),8C* ) . CC) .DC*) 

IF Bop<& THEN I=Bop 

IF Boe=5 THEN [#4 

IF Bop=!0 THEN [5 

IF Isat=! THEN 

Hs =EXP(A( I) +B( I) *LOG(X)) 

ELSE 

Hs*EXP(C( I) +D( I) .06(Х)) 

END IF 


RETURN Hs 

FNEND 

DEF FNPoly(X) 

P AC10,10),C(10).B(4) ,Nop,Iprnt,Opo,llog,Ifn,Ijsoin.Njoin 
Poly=B(0) 

FOR I=1 TO Nop 

IF Ilog*1 THEN X1-2L0GCX) 

Poly=Poly+B(I)*X1 I 

NEXT I 

IF Ilog-! THEN Poly-EXP(Poly) 

RETURN Poly 

FNEND 

SUB Poly 

DIH R(10),S(10).Sv(12),Sx(12).Xx(100).Yy(100) 

СОМ /Ср1у/ А(10.10) .С(10) .В(а4) ‚М, Іргпі, Оро, І! од, Іѓп., Іјо: п, Мјоіп 
СОМ /Ххуу/ Хр (25) , Үр(25) 

FOR I*0 TO 4 = 

ВЕТ) =9 

NEXT I 


BEEP 

INPUT "SELECT (0sFILE,!*KEYBOARD.2-PROGRAM)",Im 
Im*Im-*! 

BEER 

INPUT "ENTER NUMBER OF X-Y PAIRS",Moe 

IF Ims! THEN 

BEEP 

ШЕ "ENTER DATA FILE NAME",D file$ 

INPUT “LIKE TO EXCLUDE DATA PAIRS (1#Y,0=N)?".Ied 
ІР Іеа=1 THEN 

EE 


BEEP 

о NUMBER GF PAIRS 10° BE EXCCUDEDT IRER 
ASSIGN File TO D_files 

ELSE 


BEER 

INPUT “WANT TO CREATE A DATA FILE (1=Y,0=N)?",Yes 
IF Yes=1 THEN 

ВЕЕР 

INPUT “GIVE А NAME FOR DATA FILE",D files 

CREATE BDAT D_files.S 

ASSIGN 3File TO D_files 

END IF 

END [Р 


1-ГО 


ww— lu I 


3640 
3643 
3645 
3649 
3652 
5555 
3658 
3661 

3664 
3667 
3670 
35/3 
3676 
3679 
3682 
3685 
3688 
3631 

3634 
5597 
3700 
3703 
3706 
3703 
2212 
S715 
2718 
3721 

3724 
8727 
E730 
274216 
3736 
3739 
3742 
374S 
3748 
9751 

3754 
3/97 
3760 
3763 
3766 
3769 
2172 
4/75 
3778 
3/81 

3784 
3/87 
31130 
ЭЗ 
3736 
3739 
3802 
3805 
3808 
2911 
3814 
8817 
3820 


ВЕЕР 


INPUT "ENTER THE ORDER OF POLYNOMIAL",N 


FOR I*0 TO N«2 
Sy(I)=0 
Sx(I)=0 
NEXT I 


IF Ied-! AND Im-! THEN 


РОК [=1 ТО Ipex 
ENTER ӘҒ11е;Х,Ү 
NEXT I 

END IF 

FÜR [=] ТО Мр 
IF Im=1 THEN 


IF Opo=2 THEN ENTER @File:X,Y 
IF Opo<2 THEN ENTER SFile;Y,x 


IF Opo=1 THEN Y 
IF Iloa-! THEN 
ЈЕ Оро=2 THEN X 
Ха106(Х) 

IF Opo=2 THEN Y 
IF Opo<2 THEN Y 
END IF 


END IF 
IF Im=2 THEN 
DEEP 


2Y/X 
t=X/Y 


*LOGCXt) 
*LOG(CY) 


INPUT "ENTER NEXT X-Y PAIR",X.Y 
IF Үеѕ=1 THEN OUTPUT @File;X.Y 


END IF 

IF Im<3 THEN 
Xx (I) =X 
Yy(I)=Y 

ELSE 

Х=Хр (1-1) 
ҮзҮр(І-1) 

END IF 

R(0)=¥ 

Sy (0)=Sy(0)+Y 
S (1) =X 
Sx(1)=Sx (1) +X 
РОК Ј=1 TO N 
RCJ) RCJ- 1) *X 
Sy(J) *Sy CJ) *RCJ 
NEXT J 


FOR J=2 TO М»? 
S(J)=S(.J-1)*X 
Sx(J)=Sx(J)+S(J 
NEXT J 

NEXT I 

IF Yes-1 AND Im 
BEEP 


PRINT USING "12 
END IF 


5х(0) «Мо 

FOR I=0 TO N 
C(I) =Sy(I) 

FOR J=0 TO N 
ACT.J)=Sx(I+J) 
NEXT J 

NEXT I 

FOR I=) TO N-1 
CALL Divide(T) 


) 


) 


22 ҮНЕМ 


X,DD."" X-Y pairs uere stored in file "",108"^5;Np.D files$ 
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382 CALL Subtract(I+1) 

382 NEXT I 

3829 ВМ) См) /АСМ, М) 

3832 FOR I-s0 TO М-1 

3835 8(М-1-1)-С(М-1-1) 

3838 FOR J=0 TO I 

3841 B(N-1-I)=B(N-1-I)-R(N-1-I.N-J)*B(N-J) 

3844 NEXT J 

3847 8(М-1-1)2-8(М-1-1)/а(М-1-1,М-1-1) 

3850 NEXT I 

3853 !PRINTER IS 201 

3856 !PRINT B(»*) 

3859 !PRINTER IS 705 

3862 IF Iprnt-0 THEN 

3865 PRINT USING (eae A NEM COEFFICIENT*"" 

3868 FOR 40 ТО М 

3871 PRINT USING "15X,DD.5X.4D.7DE" ; I, BCID 

3874 NEXT I 

3877 PRINT ~ 3 

2880 PRINT USING "12X,""DATA POINT X y Y( CALCULATED) DISCR 
EPANCY"”” 

3883 FOR I=! TO Np 

3886  Yc-B(0) 

3888 ҒОК 1»! TO N 

3892 Yc=Yc+B(.J)*Xx(I) J 

3895 NEXT J 

3898 D=Yy(I)-Ye 

3901 PRINT USING "1SX,3D.4X.4(MD.SDE.,1X)":;I,Xx(I) ,Yy(I),Yc,D 

3904 NEXT I 

3907 END IF 

3910 ASSIGN ®File TO + 

3913 SUBEND 

3916 SUB Divide(M) 

3919 СОМ /Со|у/ 410,10), 27 8(4),М, Іргпе, Opo.Iloa.Ifn,IJjoin.Njoin 
3922 FOR I=" TO N 

3925 Ао=4(1І,М) 

3928 FOR J-M TO N 

3331 A(I.J)24(I,J)/Ro 

3934 NEXT J 

3997 C(I)=C(I)/ño 

3940 NEXT I 

3943 SUBEND 

3946 SUB Subtract(K) 

3949 СОМ /Coly/ 8(10,10),C(100,BCA) ,N, Iornt,Opo,I1og, Ifn.I;oin.Njoin 
3952 FOR I-K TON 

3955 FOR 14К-1 TO N 

3958 А(І.Ј) =4а(К-1,Ј) -А(І.Ј) 

3961 NEXT J 

3964 C(I)=C(K-1)-C(I) 
3967 NEXT I 

3970  SUBEND 

3973 SUB Plin 

39/6 СОМ /Cply/ A(10:10).C 
3979 COM /Xxyy/ Xx(25),Yvy( 
3982 PRINTER IS 705 
3985 BEEP 

3988 INPUT “WANT TO PLOT Uo vs Vw? (1#Y,02N)",Iuo 
3991 IF Iuo=0 THEN 

3994 BEEP 

3997 INPUT “SELECT (Q2h/h0% same +цЬе, 1 =Һ(НЕ ) /Һ(ѕт) ", Іг 


00 813 N.Iprnt Opo ss анин 
) 5 
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= ni s 


4000 
4003 
4006 
4009 
4012 
4015 
4018 
4021 
4024 
4027 
4030 
4033 
4036 
4033 
4042 
40 45 
4048 
4051 
4054 
4057 
4060 
4063 
4066 
4063 
4072 
4075 
4073 
4081 
4084 
4087 
4090 
4093 
4096 
4099 
4102 
4105 
4108 
4111 
4114 
4117 
4120 
4123 
4125 
4129 
4132 
44:35 
4138 
4141 
4144 
4147 
4150 
4153 
4156 
4159 
4162 
4165 
4168 
4171 
4174 
4177 
4180 


BEEP . 
INPUT "SELECT h/h RATIO (1*FILE.0-COMPUTED) ", Ihrat 
IF Ihrat#0 THEN 


BEEP 
INPUT “WHICH Tsat (1°6.7,02-2.2)",Isat 
END IF 


Xmin=0 

Хтах 10 

Xstep=2 

IF Irt=0 THEN 

Ymin=0 

Ymax=1.4 

Ystep=.2 

ELSE 

Ymin=0 

Ymax*15 

Ysteps5 

END IF 

ELSE 

Оро=2 

Ymins0 

Ynax*12 

Ystep=3 

Xmin=0 

Хтах 4 

Xstep*1 

END IF 

IF Ihrat-! THEN 

Ymin=0 

Ymax=18 

Ystep=3 

Xmin=0 

Хтах=9 

Хз $ер=2 š 
END IF 

BEEP 

PRINT "IN:SP1:IP 2300.2200,3300.,5800;"“ 
UM "оС 0.100.0,100: TL 2,0;" 
Sfx*100/(Xmax-Xmin) 
Sfy*100/(Ymax-Ymin) 

PRINT "PU 0,0 PD" 

FOR Xa*Xmin TO Xmax STEP Xstep 
X*(Xa-Xmin)*Sfx 

ЕНШІ ЗРАК, "05 ХТ: " 

МЕХТ Ха 

PRINT “PA 100,0:PU:” 

PRINT SPU PA 0,0 PD” 

FOR Ya=Ymin TO Ymax STEP Ystep 
Y*(Ya-Ymin) *Sfy 

АШ РВ 0. "у, "УТ" 

NEXT Ya 

PRINT “PA 0.100 TL 0 2" 

FOR Xa=Xmin TO Xmax STEP Xstep 
X=(Xa-Xmin) *Sfx 

BSNISPA':X." 100: XT" 

NEXT Xa 

PRINT “PA 100.100 PU PA 100.0 PD" 
FOR Ya*Ymin TO Ymax STEP Ystep 
Y=(Ya-Ymin) *Sfy 

ОГТ "РО Ра 100.°,Y,°YT" 

NEXT Ya 
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4183 
4136 
4183 
4192 
4195 
4198 
4201 
4204 
4207 
4210 
4213 
4215 
4219 
4222 
4225 
4228 
4231 
4234 
4237 
4240 
4243 
4246 
4249 
4252 
4255 
4258 
4261 
4264 
4267 


4270 
8273 
4276 
4279 
4282 
4235 
4288 
4231 
4234 
4237 
4300 
4303! 
4206 
4303 
4312! 
4315 
4318 
4321 
4324 
4327 
4330 
4333 


4339 
4342 
4345 
4348 
4351 
4354 
4357 


PRINT "PA 1090.100 PU” 

PRINT "РА 0. 2 56 150286 

FOR Xa=Xmin TO Xmax STEP Xstep 
X*(Xa-Xmin)*Sfx 

РЕШЕ "РА OPES" 

IF Iuo-0 THEN PRINT "СР -2.-1: 8" :Ха:““ 
IF Iuos! THEN PRINT "CP -1.5,-15;LB": Xa; "" 
NEXT Xa 

PRINT "PU PA 0.0" 

FOR Ya=Ymin TQ Ymax STEP Ystep 

IF ABS(Ya)<1.E-S THEN Ya=0 
Y=(Ya-Ymin)*Sfy 

PRINT ”РА 0,”:Ү,”” 

IF Iuo*0 THEN PRINT "CP -4,-.25;LB":Ya;"" 
IF Iuo*! THEN PRINT "CP =3,-.25: 8": та: "7" 
NEXT Ya 

XlabelS-2"OQil Percent" 

IF Iuo=0 THEN 

IF Irt=0 THEN 

YlabelS="h/h0%" 


YlabelS2"h/hsmooth" 

END IF 

PRINT "SR 1.5.2:° РА 50.-10 СР":-СЕМСХ[аве 15) /2: 70: 8": Хтаве 15 WEB 

PRINT "РА -11,50 СР 0,":-(ЕМ(Х]1ађе 15) /2+5/6: "01 0,1: 8": (1ава M 

PRINT "CP 0,0* 

ELSE 

PRINT "SP0;:SP2* 

PRINT "SR 1.2.2.4;PU PA -8,35:DI 0.1;LBU:PR 1,0.5:LBo:PR -1,0.5;LB (kH/m 


PRINT “РВ--1,0.5:58В 1.1.5:182:58 1.5.2:РК .5,.5:э .2228 .S.0: LBK) 
PRINT "РА 42,-10;0І 1,0:LBV:PR .4,-1;LBu;PR 1,.5;LBCGn/s)" 

PRINT “5РО:5Р1” š 

END IF 

Ірп=0 


EEP 

INPUT *MANT TO PLOT DATA FROM A FILE (12Y,0-2ND0?".,0kp 
Icn=0 

IF QOkp-1 THEN 

BEEP 


[INPUT "ENTER THE NAME ОЕ IHE- DATA FILES, D Fries 
IF Iuo=0 THEN 

BEEP 

INPUT "SELECT (O-LINERR. 1-LOG(X, Y)", Iloa 

END IF 

ASSIGN 9*File TO D files 


D "ENTER THE BEGINNING RUN NUMBER” ‚Md 


INPUT "ENTER THE NUMBER QF X-Y PAIRS STORED",Npairs 
IF Iuo=0 AND Ihrat-0 THEN 


BEEP 

INPUT “ENTER DESIRED HEAT FLUX" ,Q 

END IF 

BEEP 

PRINTER IS 1 

PRINT USING "4X,"*"Select a symbol:””” 

PRINT USING 14X," T ‘Star 2 Plus sigi ie 
PRINT USING "4х, "3 Circle 4" Squdre Ww» 
PRINT USING "4Х.""5 Коти 
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4360 
4363 
4366 


PRINT USING “4X.""6 
PRINT USING “4X,""7 
INPUT Sym 

PRINTER IS 705 
PRINT “PU DI" 

IF Sym=1 THEN PRINT 
IF Sym=2 THEN PRINT 
IF Sym=3 THEN PRINT 
Мп=4 

IF Iloas! THEN Nn*! 
IF Md>1 THEN 

FOR I=} TO (Md-t) 
ЕМТЕКН ФЕ / је: Ха, Та 
NEXT I 

END IF 

IF Ihrat=0 THEN 
01-0 


Rignt-side-up triangle 


10 90 #6 


Up-side-down triangl 


"G M* 
"SM+" 
"SHo" 


IF Iloa-1 THEN Q=LOG(Q) 


END IF 
FOR [=i TO Npairs 


IF Iuos0 AND Ihrat-0 


ENTER 3File;Xa,B(*) 
Ya=B(0) 

FOR Kat TO Nn 
Ya=YatB(K)*Q°K 
NEXT K 

END IF 


THEN 


IF Iuos! OR Ihrat=! THEN 


ENTER ӨР: е:Ха,Та 


IF Iuo*! ТНЕМ ҮааҮа/1000 


END IF 


IF Iuo-s0 AND Thrat=9 


THEN 


IF 110441 ТНЕМ YasEXP(Ya) 


IF Ilog-0 THEN Ya-Ql/Ya 


IF Irt=0 THEN 
IF Xa=0 THEN 
Үо-ТҮа 

Ха»! 

ELSE 

YazfYa/Yo 

END IF 

ELSE 


Hsm*FNHsmooth(Q,Xa.Isat) 


Ya=Ya/Hsm 

END IF 

END IF 

Хх(1-1)=Ха 
Yy(I-1)=Ya 
X=(Xa-Xmin)*Sfx 
Y=(Ya-Ymin) *Sfy 

IF Sym>3 THEN PRINT 
IF Sym<4 THEN PRINT 
PRINT “PA",X.Y.*" 
IF Sym>3 THEN PRINT 
IF Sym=4 THEN PRINT 
IF Sym25 THEN PRINT 
IF Sym=6 THEN PRINT 
IF Sym=7 THEN PRINT 
NEXT I 

BEEP 

ASSIGN File TO + 


то = 

"SR 1.4,2.4" 

"SR 1.2.1.6" 
о об Е. 


39. 60 
pcs 3 99529:9 5.0 


D 


= 


нов 


4543 
4546 
4549 
4552 
4555 
4558 
4561 
4564 
4567 
4570 
4573 
4576 
4579 
4582 
4585 
4588 
4591 
4594 
4597 
4600 
4603 
4606 
4609 
4612 
4615 
4618 
4621 
4624 
4627 
4630 
4633 
4636 
4639 
4642 
46 45 
4648 
4651 
4654 
4657 
4660 
4663 
4666 
4669 
4672 
4675 
4678 
4681 
4684 
4687 
4690 
4633 
4686 
4699 
4702 
4705 
4708 
3711 
4714 
4717 
4720 


END IF 
PRINT «РУ Site 
BEEP 


INPUT “WANT TO PLOT A POLYNOMIAL (1=Y,02N)?°,QOkp 
IF Okp=1 THEN 

BEEP 

PRINTER IS 1 

PRINT USING "AX.""Select line type:'""* 


PRINT USING "*8X.*"0 Solid line""*" 

PRINT USING “6X.""1 Dashed""" 

PRINT USING 76Х,772,,,5 Гопдег line - dash'""* 
INPUT Ipn 


PRINTER 15 705 


ВЕЕР 

INPUT "SELECT (QsLINEABR,1*LOG(X, Y))* , Iloa 
Iprnt=! 

CALL Poly 

IF Ічо=1 ТНЕМ 

ВЕЕР 

INPUT “DESIRE TO SET X Louer and Upper Limit (12Y,0-ND?",Ixlim 
IF Ixlim=0Q THEN 

X110 

Xul=7 

END IF 

IF Ixlim=1 THEN 


BEEP 
d "ENTER X Louer Limit",X11 


E 

INPUT "ENTER X Upper Limit",Xul 
END IF 

END IF 

FOR Xa-X11 TO Xul STEP Xstep/25 
Ісп= Ісп+1 

Ya-FNPoly(Xa) ы 
IF Iuo*! THEN Ya-Ya/1000 
Y*(Ya-Ymin)*Sfy 

X*(Xa-Xmin) *Sfx 

IF Y«0 THEN Y=9 

ЇЕ Ү»100 ТНЕМ 0010 4637 

Ри«д 

IF Ipn=! THEN Idf=Icn MOD 2 

IF Ipn=2 THEN Idf=Icn MOD 4 

IF Ipn*3 THEN Idf*Icn MOD 8 

IF Ipn*4 THEN Idf-*Icn MOD 16 

IF Ipn=5 THEN Idf=Icn MOD 32 

IF 14641 ТНЕМ Риа! 

IF Pu=0 THEN PRINT “РА. (157207 
IF Pus! THEN PRINT °PA".X,Y, “PU” 
NEXT Xa 

PRINT "РУ" 

GOTO 4285 


INPUT “WANT TO QUIT (1#Y,0=2N)?" ,ITquit 
IF Iauit*! THEN 4711 

6070 4285 

PRINT “PU SPO” 

SUBEND 

SUB Stats 

PRINTER IS 701 
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723 
4726 
4729 


4735 
4738 
4741 
4744 
4747 
4750 
4753 
4756 
4753 
4762 
4765 
4768 
4771 
4774 
4777 
4780 
4783 
4786 
4789 
4732 
4795 
4798 
4801 
4804 
4807 
4810 
4813 
4816 
4819 
4822 
4825 
4828 
4831 
4834 
4837 


4840 
4343 


З з) 

К =!) 

ВЕЕР 

IF Iplot=! THEN ASSIGN 8File ТОР Ғі1е% 

BEEP 

INPUT "LAST RUN No?(0=QUIT)” Nn 

IF М№п=0 THEN 4849 

Nn=Nn-J 

Sx=0 

Sy=0 

51-0 

Sxs=0 

Sys=0 

Szs=0 

FOR I=! TO Ма 

J=J+1 

ENTER @File:Q.7T 

H=Q/T 

5х20х-0 

Sxs=5xs+Q°2 

Sy*Sy+T 

Sys=Sys+T 2 

Sz=Sz+H 

525=5;5+Н'2 

NEXT I 

Qave=Sx/Nn 

TavesSy/Nn 

Have=Sz/Nn 

Sdeva=SQR( ABS ( (Nn*Sxs-Sx°2)/(Nn*(Nn-1)))) 

Sdev t» SQORCABS((Nn*Sys-Sy 2) /CNn*(GNAn-122)) 

Sdevh* SORCABSC(ONn*Szs-Sz'2)/(Nn*(iNn- 12220) 

Sh=100*Sdevh/Have 

Sq=100*Sdevq/Qave 

St=100*Sdevt/Tave 

IF K=1 THEN 4843 

PRINT 

AINT USING "11X,""DATA FILE:"",!4R";FileS | 
N 

PRINT USING "11X,""RUN Htube Саен Qdp SdevQ Thetab SdevT""* 


Kal | 
PRINT USING "11X,DD.2(2X,D.3DE. 1X,3D.2D) ,2X,DD. 3D, 1X.3D.2D* ; J, Have, Sh, Qave 


Поа. Таме .5: 


4846 
4849 
4852 
4855 
4858 
4861 
4864 
4867 
4870 
4873 
4876 
4879 
4882 
4885 
4888 
4831 
4834 
4897 


0010 4735 

ASSIGN #Filei TO >» 

PRINTER IS 1 

SUBEND 

SUB Coef 

ШЫ РУ 8(10,10) ,СС10),8(4) ,М, ХргпЕ Оро, (109, [Ееп, Х јо1п,МјоЈп 
INPUT “GIVE A NAME FOR CROSS-PLOT FILE",CofS 
CREATE BDAT Cef$.2 

ASSIGN @File TO Cpfts 

BEEP 

заці "SELECT (O=LINEAR.1=LOG(X.Y))"*,Ilag 
INPUT "ENTER OIL PERCENT (-1=STOP)",Boo 

ТЕ Вор<0 ТНЕМ 4900 

CALL Poly 

OUTPUT @File;Bop.B(*) 

GOTO 4882 


m 


4900 
4903 
4906 
4909 
4912 
49151! 
4918 
4921 
4924 
4927 
4930 
4933 
4936 
4939 
4942 
4945 
4348 
4951 
4954 
4957 
4960 
4963 
4966 
4969 
4972 
4975 
4978 
4981 
4984 
4987 
4990 
4993 
4996 
4999 


5005 
5008 
5011 
5014 
5017 
5020 
5023 
5026 
5029 
5032! 
5035 
5038 
5041 
5044 
5047 
5050 
5053 
5056 
3059 
5062 
5065 
5068 
5071! 
5074! 
5077 


ASSIGN sFile TO * 

SUBEND 

SUB Hilson(Cf,Ci) 

COM ЖАШ реси Во ПРИ МЕН 

DIM Emf (12) 

EGE, PLOT SUBROUTINE DETERMINE CF AND CI 
INPUT “ENTER DATA FILE NAME" ,Files 

BEEP 


PRINTER IS 1 

PRINT USING "4X.""Select option: '"" 
PRINT USING ”"4X.”= -0 Vary CF and CiT” 
PRINT USING "4X."" | Fix Cf Vary Ci""" 
PRINT USING *"AX."" 2 Vary Cf Fix Ci""" 
INPUT "ENTER OPTION",Icfix 

PRINTER IS 7201 

IF Ictix=0 THEN 4960 

IF Icfix>0O THEN BEEP 

IF Icfix=! THEN INPUT "ENTER Cf" .Cst 
IF Icfix-2 THEN INPUT “ENTER CI”,Ci 
PRINTER 15 1 

INPUT “Want To Vary Coetf?(1=Y,02N)",Iccoef 
IF Iccoef=1 THEN INPUT "ENTER COEFF”,R 
PRINTER IS 701 

IF Icfix=0 OR Icfix=2 THEN Cfa=.004 

IF Icfix=1 THEN Cfa=Csf 

Cia=Ci 

Xn=.3 

Fr=.3 

4140 

Rr=3. 

IF Iccoet=! THEN Rr=R 

PRINTER IS 1 

PRINT Do.0i.kKeu 

ASSIGN ®File TO Files 

ENTER sFile:Nrun,Dates.Ldtcl,Ldtc2,Itt 
Rus0o*LO0GCDo/D1)/(2*Kcu) 

Sx=9 

Sy 20 

Sxy 29 

5х2з0 

Sy2=0 

FOR [=1 ТО Мгип 

ENTER 93File:Bop.ToldS.Emt(*) Fms 
CONVERT EMF'S TQ TEMPERATURE 

РОВ Je! 70 5 

TCJ) *FNTvsv CEmf CJ)) 

NEXT J 

Твасч(Т(1)-1(2))».5 

Тама- Ї (5) 

Grad-37.9853*.104388*Tava 
Tdropsztmt(7)*1.E*6/(10.*Grad) 
Тамас» Т(5) -Тагоре.5 

IF ABS(Tava-Tavac)».01 THEN 
Tava*(Tava*Tavqc) *.5 

GOTO 5050 

END IF 


Compute properties of water 
КшчҒМКш( Тама) 
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9093 
5236 
5233 
5242 
5245 
5248 
5251 
5254! 
5257 


МишааЕММиш( Тамач) 
Cew=FNCow( Tavg) 
Ргш=Е МР ги ( Тауа) 
Rhouy*FNRhou(Tavqa) 


1 Compute properties of Freon-l14 


Lmatd=Tdrop/LOG((T(5)-Tsat)/(1T(S)-Tdrop-Tsat)) 
IF Jj*0 THEN 
Tu*Tsat*Fr*Umntd 
Thetab=Tu-Tsat 
3141 

ENDO IF 
Tf-*(CTu*Tsat)*.5 
Rho=FNRho( Tf) 
Mu=FNMu( Tf) 
K=FNK(TF) 
Cp=FNCp( TF) 
Beta=FNBeta(Tf) 
Hf g=FNHfg(Tsat) 
Ni=Mu/Rho 
Alpha-K/(Rho*Cp) 
Pr=Ni/Alpha 


Analysis 

COMPUTE MDOT 

A-PI*(Do'2-0i1'2)/4 

Р=РІ+0о 
Mdot-23.9657E-3*Fms*(3.61955E-3-Fms*(8.32006E-6-Fms*(1.23688E-7-Fms*4.31897 
) 


Q-Mdot*Cpu*Tdrop 

COMPUTE NATURRL-CONVECTIVE HEAT-TRANSFER COEFFICIENT 
FOR UNENHANCED END(S) 

НБаг= 190 

Fe*(Hbar*P/(Kcu*0))'.5*Lu ” 

Tanh=FNTanh(Fe) 

Theta*Thetab*Tanh/Fe 

Xx=(9.81*Beta*tThetab*Do’ 3*Tanh/(Fe*Ni*Alpha)) °.166667 
Уу=(1+(.559/Рг; * (9/16) ) * (8/27) 

Hbarc*K/Do*(.6-.387 *Xx/Yy)^2 

IF ABS( (Hbar-Hbarc) /Hbar)>.001 THEN 
Hbar=(HbartHbarc)+*.5 

GOTO 5176 

END IF 


! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED ENDS 


Q1-2(bar*P*Kcu*2) ' .5*«Thetab*Tanh 

Qcz20-2*01 

Rs=PI*D2*L 

COMPUTE ACTUAL HEAT FLUX 

Qdp=0c/As 

IF Ictix=0 OR Icfix>1 THEN Csf=1/Cf*(1t./Rr) 
Thetab=Csf /Cpo*#Hfg*(Qdp/ (MutHfq) *(.014/(9.381*Rho))°.5)°(t/Rr) *Prot.7 
Ho=Qdp/Thetab 

Omeqa=Ho/Cf 

Оо / (Р Т»роч, 4 під) 

VysMdot/(Rhou*PI*01'2/4) 

ReyszRhouw*Vu*b i/Muua 

TuizTurQ*Ru/(PI*0os*t) 
Gama-Ku/Di*Reu^.9*Pru' (1/3.) *(Muua/FNMuut(tTu1)) ^ .14 
PRINTER IS 1 

Yu=(1./Uo-Rw) *Omega 
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5437 


Хшедпеааєдо/ (батає0 1) 

Sx=Sx+Xw 

Sy *Sy*Yu 

Sxy=Sxy+*Yu*Xu . 
Sx2=Sx2+Xu*Xu Миг 
Sy2=Sy2+Yu*Yu 

NEXT I 

ASSIGN SF. ile TO + 
M=(Sx*Sy-Nrun*Sxy)/(Sx*Sx-Nrun*Sx2) 
C»(Sy-Sx*M) /Nrun 

IF Icfix*0 OR Icfix=3 OR Icfix=4 THEN 
Сіс=1/М 

Сїса1/С 

END IF 

IF Icfix-1 THEN 

С1с=1/М 

Cfc=Cf 

END IF 

IF Icfix=2 THEN 

Cie=Ci 

Сїс=1/С 

END IF 

IF ABS((Ci-Cic)/Cic)>.001 OR ABS((Cf-Cfc)/Cfe)>.001 THEN 
Сіз(СізСіс)ж.5 

Cf*(Cf*Cfc)*.5 

PRINTER IS 1 

PRINT USING "2X."" Csf * "".MZ.3DE.2X."" Ci = "*,MZ.SDE'-CS OMNEM 
PRINTER IS 701 


GOTO 5002 

END IF 

PRINT 

PRINTER IS 701 

PRINT USING 723X." Cf UID 

PRINT USING "3X,""ASSUMED "ЛА ЗВЕ ЗА НАЗОВЕ G I a Ы 


ERINI USING "3X." “CALCULATED "" ,Н2.ЗрЕ.ЗАХ,И222 0Е70-122271 
N 

Sumz-*Sy2-2*M*Sxy-2*C*Sy *M^2*Sx2*2*M*C*Sx*Nrun*C^2 
PRINT USING “10Х, ""бит оф Squares = "^,Z.3DE":Sum2 
PRINT USING "10X.""Coetficient » "",D.DDD":Rr 

SUBEND 

DEF ЕММиш( Т) 

A2247.9/(T*133. 15) 

Muz2.4E-5*10'A 

RETURN Mu 

FNEND 

DEF FNCoutT) 
Cpu24.21120858-T*(2.26826E-3-T*(4.42261E-5*2.71428E-7 *T)) 
RETURN Cpu*1000 

FNEND 

DEF FNRhou(T) 

Ro-999.52946*T*(.01269-T*(5.482513bE-3- T*1.234147E-53)) 
RETURN Ro 

FNEND 

DEF FNPru(T) 

Pruz*NCou( T) *FNMuu(CT) /;FNKu(CT) 

RETURN Pru 

FNEND 

DEF FNKw(T) 

Х=(1+273.15)/273.15 
Кш=-,92247+Х*(2.3395-Х*(1.3007-Х*(.52577-.07344*Х))) 
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5440 
5443 
5446 
5449 
5452 
5455 
5453 
БАИ 
5464 
5467 
5470 
5473 
_ 5476 

5479 
5482 
5485 
5483 
5491 
5494 
5447 
5500 
5503 
5506 
5509 
5512 
5515 
5518 
5521 
5524 
5527 
5530 


RETURN Kw 

FNEND 

SUB Plot 

COM /Cply/ 8A(10,10),Ct 100,B( A) ,Nop, IDrnt,O0po.Ilog.Ifn,Ijoiíin.Njoin 
DIM Bs(3) 

INTEGER Ii 

PRINTER IS 1 

BEEP 

P "SELECT HEATING MODE (O=ELECTRIC, t=WATER)*,Ihm 
dv= 

BEEP 

INPUT "LIKE DEFAULT VALUES FOR PLOT (1#Y,02N)?°,Idv 


PRINT USING “4X,"*"Select Орбіоп:""" 
PRINT USING °6X.""0 q versus deita-T""" 
PRINT USING "6X.""1 h versus deita-T""" 
PRINT USING *65X,""2 h versus q"''" 

[МРОТ Оро 


INPUT "SELECT UNITS (Q-»SI,17ENGLISH)",Iun 
PRINTER IS 705 

ТЕ Idv<>1 THEN 

BEEP 


ЕНТ REDITERONUMBERCOF. CYCLES FOR X-AXIS",Cx 
BIBSICCENTERSNUMBER OF CYCLES РОК 7-ах15" ,Су 


INPUT "ENTER MIN X-VALUE (MULTIPLE OF 10)".Xmin 
Mime ENTER HIN -VALUE (MULTIPLE OF 10)".Ymin 


IF Opo=0 THEN 
IF Ihm=1 THEN Б B 
Су з | 

Cx=2 

Xmin= | 
Ymin=10000 
ELSE 

Cy 23 

Сх=3 

Хтап=. | 

Утуп= 100 

END IF 

END IF 

IF Opo=1 THEN 
IF Ihm=! THEN 
Су=2 

Сх=2 

Xmin=.1 
Ymin=100 

ELSE 

Cy=3 

Сх=2 

Хпіп” 1000 
Хтіпе100 

END IF 

END IF 

IF Opo=2 THEN 
IF Ihm=1 THEN 
Cy=2 


2253 


Cr» 

Xmin=10000 

Ymin=100 

ELSE 

Су «З 

Сха3 

Xmin*.1 

Ymin=100 

END IF 

END IF 

END IF 

BEEP 

PRINT "IN;SP1;IP 2300,2200.8300,6800;" 
PRINIOSSG 0:200.0, 1004 TL 727025 
Sfx=100/Cx 

Sfy2100/Cy 

BEEP 


INPUT “WANT TO BY-PASS CAGE? (12Y,0-N)",Ibyp 
IF Ibyp=! THEN 5049 

PRINT "РУ 0,0 РО" 

Ма» З 

FOR I*1 TO Cx! 

Хаё= Хм: п+10* (1-1) 

IF I2Cx*1!1 THEN Nn=1 

FOR Js! TO Nn 

IF J*1 THEN PRINT "TL 2 0" 
[Е J=2 THEN PRINT "TL 1 0" 
Xa=Xat*J 
X*LGT(Xa/Xmin) *Sfx 

PRINT "РАМЕ ОУ Т" 
МЕХТ Ј 

NEXT I 

PRINT "РА 100,0:Р09: 

PRINT "PU PA 0,0 РО" Е 
Nn*9 ` 

FOR I2»! TO Cy-*1 
Yat*Ymin*! O^ (CI-1) 

IF I=Cy+! THEN Мпз1 

FOR J=1 TO Nn 

IF J=1 THEN PRINT “TL 2 0" 
IF J=2 THEN PRINT "TL 1 O° 
Үа=Үзё+Ј 
Y=_GT(Ya/Ymin) *Sfy 

PRINT РА О Y. YT" 

МЕХТ Ј 

МЕХТ 1 

РВ "PA 0:000 IL 0 25 
Мп=9 

FOR І=1 ТО Сх+1 
Xat*Xmin*10^ CI-1) 

IF ЇаСх-! ТНЕМ Маа! 

FOR J=1 TO Nn 

IF J=1 THEN PRINT *"TL 0 2" 
IF J>1 THEN PRINT "TL 0 1" 
Xa=Xat*J 
X*LGT(Xa/Xmin) *Sfx 

PRINT А РА C DOS IS 
NEXT J 

NEXT I 

PRINT "PA 100.100 PU PA 100,0 PD" 
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5803 
5806 
5809 


Мп=9 

FOR [=] ТО Су! 

Yat=Yminel0°(I[-1) 

IF [=Cy+1!1 THEN Nn*! 

FOR J=1 TO Nn 

IF J=! THEN PRINT “TL 0 2" 

IF J>1 THEN PRINT “TL 0 1" 

Ya=Yat*J 

Y=LGT(Ya/Ymin) *Sfy 

PRINT "PD PA 100,",Y,"YT" 

NEXT J 

NEXT I 

PRINT “PA 100.100 PU" 

PRINT "РА 0.-2 55 1,5,2” 

IizLGT(Xmin) 

РОК 1=1 ТО Сх+1 

Xa=Xmin*10 (I-I) 

X=LGT(Xa/Xmin)*Sfx 

PRINT РАТ. "0: “ 

PERSO THEN PRINT “СР -2,-2:1810:РЕ -2.2: 8" :11::“" 
ПЕКИ ОО НЕН РЕ СР -2,-2:1В10:РЕЁ 0,268" 11: 
18538952 

МЕХТ 1 

PRINT "PU PA 0.0" 

Ti=LGT(Ymin) 

Y10=10 

FOR =1 ТО Су+1 

Ya=Ymine10° (I-17) 

Y= nM o M ty 

PRINT "РА 0, 

PRINT "СР -4,- 225: FD IO PRPO O OD DB" 01 
Гізі1%1 

NEXT I 

BEEP I 
INPUT "MANT USE DEFAULT LABELS (1=Y.,0=N)?".Idl 
IF Idl<>1 THEN 

BEEP 

EJT "ENTER X-LABEL",XlabelS 


INPUT "ENTER Y-LABEL",YlabelS 

END IF 

IF Opo<2 THEN 

PRINT "SR 1,2:?U PA 40,-14:" 

PRINT "LB(CT;PR -t. 07 5 PD PR 1.2,0 PU:PR .5,-4;LBuo;PR .5,1;" 
PRINT "LB-T:PR .5,-1;LBsat:PR „5,1:° 
ТЕ Iun=0 THEN 

DENIM EB) К)" 

ELSE 

PRINT 218) (F)* 

END IF 

END IF 

IF Opo=2 THEN 

IF Iun=0 THEN 


no. Е РУ P0130 3; Ba (H/m SR 1 S:IPR 0.5.1:LB2:SR 1.5.2;PR 
ELSE 

Бок „52 РО Ра 34.-14: Ва (Btu/hr:PR .5..5:LB. РВ .5,- 5547 

u ss (pan .5.1:5К 1,1.5:LB2;SR 1.5,2:PR ,S.-1:LB); 

END IF 

END IF 


IF Opo=)) THEN 


025 


IF Iun-0 THEN 
PRINT "SR 1.5.2:PU PA -12.40:DI 0,!: Ва (W/m;PR -1,0.5 5 КИКИ И 


"РА. В)” 


ELSE 

PRINT “SR 1.5.2:РУ РА - BS 32:DI 0.1:LBa (Btu/hr:PR =.5..5: 8. РК СЫ 
PRINT "LB8ft; SR. ura S PRE 1,.5: LB2; PR 1, S:SR (5 2018 :15 

END IF 

END IF 

IF Opo>0 THEN 

IF Iun=0 THEN 


|! PRINT "SR 1.5,2;PU PA -12,38:DI 0,1:L8h (M/m;PR -1,.5;SR 1,1.5;LB2;SR 1. 


25% 26 
PRINT "SR 1.2,2.4;PU PA -12,37:DI 0,1 8ћ:РК 1,025 Во: РК 1n SRL 


PRINT "PR -1,..5;SR 1.1.5;LB2;SR 1.5,2:PR .5.. 9 ЕВ. РКО eee 


ELSE 

PRINT “SR 1.5.2:PU PA -12.28:D0I O.13LBh (Btu/hr:PR -.S..S:LB.;PR 3500 
РКІМТ “ВЇЬСРБ -1,.5:58 1.1.5:182:58 1.5,2:РВ .5,.5: 2 РНГ 115. LBF)* 
END IF | 

END IF 

IF Idl*0 THEN 

PRINT "SR 1.5.2:PU PA 50.-16 CP'":-LENCXIabelS)/2: "OQ: EB" лаве! ae 
PRINT "PA -14.50 CP 0,"*;-LEN(YL1abelS)/2*5/8;"DI 0,1;1,83";YlabelS;"" 
РЕТМТ”СР” 12070129 

END IF 

Ipn=0 

Х11=+1.Е+6 

Хці=-1.Е+6 

Ісп=0 

Іп» 0 

Tyoin=} 

BEEP 

INPUT “WANT TO PLOT DATA FROM A FILE Roper. 0=N)?" ,Ok 

IF Ok=! THEN 

BEEP 

INPUT "ENTER THE NAME OF THE DATA FILE",D files 

ASSIGN 9*File TO D_files 

BEEP 

BEEP 

Do "ENTER THE BEGINNING RUN NUMBER” .Md 


СЕР 
INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED",Npairs 


I BEEP 


INPUT “CONNECT DATA WITH LINE (12Y,02N)?", [cl 
BEEP 

PRINTER IS 1 

PRINT USING "4X."*"Select a symool:""" 

PRINT USING “6Х.""1 Star 2 Plus зіап""" 
PRINT USING="6A.°°S Circle 4  Зацаге""" 
PRINT USING "5Х.“"5 Ronmous'"" 

PRINT USING "6X.“°"6 Right-side-up triangle""" 
PRINT USING °6X,°°7 Up-side-doun triangie""" 
INPUT Sym 

PRINTER IS 705 

PRINT "PU DI" 

IF Sym=1 THEN PRINT "SM*" 

IF Sym=2 THEN PRINT "5М+“ 

IF Sym=3 THEN PRINT "5Мо" 

IF Md»! THEN 

РОК І»! ТО (Ма- 1) 
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ENTER eFi1le:Ya.Xa 

NEXT I 

END IF 

FOR I=! TO Npairs 
ENTER @File:Ya.iia 

IF I=? THEN Q1=Ya 

IF I*Npairs THEN Q2=Ya 
IF Opo~t THEN Ya=Ya/Xa 
IF 0po-2 THEN 


E IF 

IF Xa«X11 THEN X1l*Xa 

IF Xa>Xul THEN Xul#Xa 

IF Iun=1 THEN 

IF Opo<2 THEN Xa~Xa*#1.8 

IF 0po»0 THEN Үа=Үа+, 1761 

IF Opo~0 THEN YazYax.317 

IF 0po-2 THEN Xa-Xa*.317 

END IF 

X*LGT(Xa/Xmin)*Sfx 

YeLGT(Ya/Ymin)*Sfy 

Кј=0 

CALL Symb(X. Y, Sym, Icl,K;) 

GOTO 6268 

IF Sym»3 THEN PRINT "SM" 

IF Sym<4 THEN PRINT “SR 1.4,2.4" 

IF Ici-0 THEN | 

PRINT “PA",X,Y,°° 

ELSE 

PRINT “РА”,Х,Ү,“РО” 

END IF | - 
IF Sym>3 THEN PRINT "SR 1.2.1.6" 
IF Sym=4 THEN PRINT 5 
IF Sym=S THEN PRINT “UC3.9.99. 
IF Sym=6 THEN PRINT "UC0.5.3.9 
IF Sym=7 THEN PRINT "0С0,-5.3, 
NEXT I 

PRINT "Р" 


BEEP 

INPUT “WANT TO LABEL? (1#Y,0"N)",Ilab 
IF Ilab-i THEN 

BRINT *SPO:SP2* 


BEEP 

IF Xlab=0 THEN 

Xlab=5 

Ylab=85 

ШРШ РЕМТЕК INITIAL X.Y LOCATIONS” ,Xlab.Ylab 
Xtt=Xlab-5 

YttzYlab*8 

PRINT *SR 1,1.S" 

BEEP 

INPUT “SELECT LABEL TYPE (O=FILE NAME, 1=TUBE TYPE)",Ifntt 
IF Ifntt=0 THEN 


6 
= 
x 
9, on 
gu 


ESI L SM:PAT tt, Ytt, "LB Z Heat File” 
ШЕС-ТІЕЗ-3 

Пном РА“ ХЕС, 055. "8 4^ Heat File" 
BLSE 

ДЕРИ "БИРА" БЕБЕ." В Z Heat Tube" 


T27 


6391 
6334 
5397 
6400 
6403 
6406 
6409 
6412 
6415 
6418 
6421 
6424 
6427 
6430 
6433 
6436 
5439 


6442! 


6445 
6448 
6451 
6454 
5457 
5460 
5463 
5466 
6463 
6472 
6475 
6473 
6481 
6484 
6487 
6490 
6493 
6496 
6499 
6502 
5505 
6508 
5511 
6514 


Еру ЕЕ 5 

Ро РЯ СА. ЕЕ. eH 0:1 Flux Type" 
END IF 

IF Sym=t THEN PRINT "SHM«" 

IF Sym=2 THEN PRINT “SM+" 

IF Sym=3 THEN PRINT “SMo" 
Klab=1 

END IF 

Kj! 

CALL Symo(Xlab.Ylab.Sym,Icl,Kj) 
PRINT "SR 1.1.5:SM^ 

ІҒ Әуп<4 THEN PRINT "PR 2,0" 
INPUT “ENTER BOP",Bop 

IF Ifntt=1 THEN 


BEEP 

INPUT “ENTER TUBE TYPE".TtypeS 

END IF 

IF Bop<!10 THEN PRINT “PR 2,0:LB":Bop;"" 
IF Bop»9 THEN PRINT "PR .5,05;1,8";Bop;"" 
[hf =0 

IF Q1>02 THEN Ihf=! 

IF Ihf=0 THEN PRINT “PR 4,0:LBInc” 

IF Ihf=t THEN PRINT “PR 4,0:L8Dec" 

IF Ifntt-0 THEN 

PRINT 2РН:22:01:35:021103:27 


ELSE 

РКИМЕЗРК 2 0: 287 Гіурез:;”” 
END IF 

PRINS SPO SPI SR 1.5.2; 
Ylab=Ylab-5 

END IF 

BEEP 

ASSIGN ıle TO e+ 
Х11еХ11/1.2 

KulsXul*1.2 

GOTO 3040 

END IF 

PRINT "PU SM" 

BEE 


p 
INPUT “WANT ТО PLOT A POLYNOMIAL (1=Y,0=N)?",Go_on 
IF Go_on=1 THEN 

BEEP 

PRINTER 15 1 

PRINT USING "4X.""Select line type:'"" 

PRINT USING "*5X."*0 Solid line""" 

PRINT USING "SX,""1 Dashed""" 

PRINT USING "5X,"'2,,,5 Longer line - dash ~~ 
INPUT Ipn 

PRINTER IS 705 


BEEP 

INPUT “SELECT (OQ=LIN.1#LOG(X.Y))"*,Ilog 
Iprnt=! 

CALL Poly 

IF Ifn=0 THEN 

Nyoin=! 

BEEP 

INPUT “ENTER NUMBER OF FILES TO JOIN?",Nyoin 
END IF 

Тјо1п=0 

IF Ifn<Nyoin THEN Ijoin-! 
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IF Ifn-0 OR Ijoins! THEN 
FOR Ij2s0 TO 3 
Bs(Ij)28s(Ij)*B( Ij) 
NEXT Ij 

Ifn=lf nel 

END IF 

IF Njoin=Ifn THEN 
FOR I; j=0 TO 3 

Bi Iy)=Bs(Iy)/Nyoin 
Ва([ј)=0 

МЕХТ 1] 

Ifn=0 

ELSE 

GOTO 6067 

END IF 


BEEP 

INPUT "ENTER Y LOWER AND UPPER LIMITS",Yll.Yul 
РОК Хх=0 ТО Сх STEP Cx/200 

Xa=Xmin*10 Xx 

IF Xa«X11 OR Xa»Xul THEN 6555 

IcnszIcn-*! 

Pu=0 

IF Ipn=! THEN Idf=Icn MOD 2 

IF Ipn=2 THEN Idf#Icn MOD 4 

IF Ipn=3 THEN Idf=Icn MOD 8 

IF Ipn#=4 THEN Idf#Icn MOD 16 

IF Ipn=S THEN Idf=Icn MOD 28 

IF Idf=1 THEN Pu=1 

IF Opo=0 THEN Ya=FNPoly(Xa) 

IF 0po-*2 AND Ilog=0 THEN Ya=Xa/FNPoly( Xa) 
IF Opo=2 AND Ilog=! THEN Ya=FNPoly(Xa) 

IF Ороз1! THEN Ya=FNPoly( Xa) 
IF Ya<Ymin THEN 6655 

IF Ya<Yll OR Ya>Yul THEN 6655 
IF Iun=! THEN 

IF Оро<2 THEN Xa-Xa*1.8 

IF 0po^»0 THEN Yas-Yas.1761! 

[Е 02о=0 THEN Ya-Ya*.317 

IF Opo#2 THEN ХазХа».317 

END IF 

Y*LGT(Ya/Ymin)*Sfy 
X*LGT(Xa/Xmin) *Sfx 

IF Y<0 THEN Y=9 

IF Y>100 THEN GOTO 6655 

IF Pu=0 THEN PRINT "РЯВ",Х.ї,"РО" 
IF Pu*1 THEN PRINT “РА”,Х, у, "РУ" 
NEXT Xx 

PRINT "PU* 

GOTO 6067 


INPUT “WANT TO PLOT REILLY^S SMOOTH-TUBE DATA? (1-7,0-М)“,Ігіу 
IF Opo=0 OR Opo=1 THEN 

Х11=3 

Хи І «20 

END IF 

IF 0po-2 THEN 

Х11=10000 

Хи!=100000 

END IF 


T29 


IF Irly=1 THEN 

ТЇ 220 

Уи | =70 

ВЕЕР : . 
INPUT "ENTER LOHER AND UPPER Y-.IMITS FOR PLOTTING",Yll.Yul 
FOR Xx=0 TO Cx STEP Cx/200 

Xa=Xmin*10 Xx 

IF Xa<X11 OR Xa>Xul THEN 6757 

X1=L0O0G(Xa) 


ЈЕ Оро=0 THEN 71=-2.54028372-1+Х1+(4,.9720151-Х 12.515 а ВИ ИЕ 
ТЕ Орогі THEN Y1-2-2.5402837E-1*X1*(3.9720151-X1*2.5124787E - 1) 
IF Qpo=2 THEN Y12-3.7073801E-1*X1*(98.7259190E-1-X1*6.3825842E 739 


Үа=ЕХР(Ү1) 

Y*LGT (Ya/Ymin)*Sfy 
X=LGT(Xa/Xmin) *Sf x 

Три=0 

IF Y«Y1l THEN Ipu=1 

IF Y>Yul THEN GOTO 6757 

IF Ipu=9 THEN PRINT *PA",X,Y, "PD" 
IF Ipu*1 THEN PRINT “PA",X.Y, *PU" 
NEXT Xx 

PRINT "РУ" 

END IF 

BEEP 

INPUT “WANT TO PLOT ROHSENORWK CORRELATION? (1=Y,0=#N)°,Irohs 
IF Irohs=t THEN 

Ylls15 


INPUT "ENTER Tsat (Deg C)*,Tsat 
Csf=.0040 

BEEP 

INPUT “ENTER Csf (DEF=0.004)"*.Csf 
Tf=Tsat+2 i 
FOR Xx=0 TO Cx STEP Cx/200 
Xa=Xmin*10 Xx 

IF Xa«X11 OR Xa»Xul1 THEN 6889 
Xl=LOG( Xa) 

IF Qpo<2 THEN Tf=TsatrXa/2 
RhozFNRho( TF) 

К=ЕМК (ТЕ) 

Mu=FNMu( TF) 

СреРМСо (ТЕ) 

НҒачҒМНҒа(Тса%) 

Ni=Mu/Rho 

Pr=Cp#Mu/K 

Omeaqa=Cst #Hfaq/Co*((.014/(9.31*Rho)) *.S/(Mu*Hfg))°(1./3) *Pr°t.7 
IF Opo=9 THEN Ya=(Xa/Omega)? 3 

IF Opo=! THEN Уа=(Ха/Отеза) `З/Ха 
IF Opo=2 THEN Үа=Ха` (2. /3) /Отеда 
IF Оро-2 ТНЕМ 

Tfez1sat*Xa/Yas.5 

IF ABS(Tf-Tfe)>.01 THEN 
Tf=(Tf+Tfc)*.S 

GOTO 6814 

END IF 

END IF 

Үзі ОТ(Үа/Үтіп) «5Ғу 
X*LGT(Xa/Xmin)*Sfx 

Ipu=0 
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7162 
7165 
7168 
7171 
7174 
227 


РВЇМТ 

PRINT USING "10X."" Fms Tin Tev Чу `2 Патар > 

ІР Кай THEN 

PRINTER IS 701 

PRINT 

PRINT " Month. date and time: ":Dates$ 

IF Itt-0 THEN PRINT USING "10X.""Tube Type: Wieland Smooth"*"" 
IF Itt=! THEN PRINT USING "10X.""Tube Type: High Flux"*" 

IF Itt=2 THEN PRINT USING "10X,""Tube Type: Turbo В” 


PRINT 
PRINT USING "*10X,"" Fms Tin Tev Уш^2 Тагор""" 
PRINTER IS 1 
Ка! 
ЕМО ЇЕ 
BEEP 
INPUT "ENTER FLOKMETER МАН Fms 
OUTPUT 709;"AR AFO AL4 V 
FOR L=0 TO 4 
OUTPUT 709;"AS Sñ“ 
х L»0 AND L«4 THEN 7141 
=0 
FOR I=0 10 9 
ENTER 709:E£ 
S*S«E 
NEXT I 
IF L=0 THEN Emt(0)=ABS(S/10) 
IF L=4 THEN Emt(1)2=ABS(S/10) 
NEXT L 
OUTPUT 709:“АВ AF20 А120 ЧНІ" 
OUTPUT 709;"AS SA” 
Etp=0 
FOR I=0 TO 9 
ENTER 709:Et - 
Etp-*Etp*Et 
NEXT I 
Etp=Etp/10 
Tin=FNTvsv(Emt (1) ) 
Tev2FNIvev(Emt (0) ) 
Grad=37.9853+.104388+*Tin 
Mdot=3.9657E-3+Fms4+(3.619SSE-3-Fms*(8.32006E-5-Fms#(1.23688E-7-Fms*4. 31897 


ESTO))) 


7180 
7183 
7186 
7183 
2132 
7135 
7138 
7201 
7204 


Vy*Mdot/(1000*PI*Di'2)*4 
Tdropz£tp*1.E*6/(10*Grad) 
БО! USING 7310Х.3(00.00.4Х),1Х.2.00.4Х,М2.40" :Епе,Тіа,Тем.Чш'2,Тагор 
КТ "WANT TO ACCEPT THIS DATA SET? (1=¥,0=N)",Ok 
= J+] 
IF Ok=0 THEN 
J=J-1 
GOTO 7102 
ЕДЕЗЕ 
QUTPUT @File:Fms.Emt(*) ,Etp 
PRINTER IS 701 
PRINT USING "10Х,3(00.00,4Х), Х,2.00,4Х.М2.40":Ёпе,Тїп.Ївм,Уш 2,Тагор 
PRINTER IS 1 
BEEP 
INPUT “WILL THERE BE ANOTHER DATA SET? (127,0-М)",бо оп 


IF бо оп”! THEN 7102 
END IF 


ASSIGN File TO + 
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гааг 
7240 
7243 
iles 
7246 
7249 
1252 
7255 
1458 
7261 
7264 
(257 
7270 
(213 
7276 
2203 
2282 
7285 
7288 
229 
7234 
7227 
7300 
7303 
2306 
7308 
2312 
7315 
7318 
732] 
7324 
2222 
7330 


PRINTER 1597 


РЕТМТ 
PRINT 


POIN ER ЗЕ 
SUBEND 


USING 


SUB Uoprt 


PRINTER IS 1 
BEEP 


"10X.""NOTE: "".77,"" data sets are stored In пен 


INPUT “Enter Uo File Name“.FileS 
BEEP 


INPUT "Number of Data Runs',Nrun 
INPUT “Do You Want a Plot File?(1=Y,0=N)" ,Iplot 


ВЕЕР 


IF Iplot*! THEN 


INPUT *Give Plot File Name",P file$ 


CREATE BDAT P file$.4 


ASSIGN @Plot TO P_files 


END IF 
PRINTER IS 701 
PRINT 


PRINT 
PRINT 


ASSIGN *File TO FileS$ 


USING 


Water Vel Чо. > 


IF Iplot-! THEN ASSIGN GFile! TO P files 


РОК 1=1 


и Те ове 
PRINT USING "15X,D.0D.6X.4Z.3DE" :Vu.Uo 


NEXT I 

ASSIGN ФЕ: іе ТО + 

ASSIGN Е: еї Т0 + 

"10Х,""МОТЕ: "7",21,"" data sets are stored іп Біїе "7 БА. И 


PRINT 


n,Filed 
IF Iplot=!{ 

7336 PRINT USING 
n.P_files 


(839 


2333 
7342 
7345 
7348 
235! 


ЕШТЕ 
РЕТТЕН IS i 
SUBEND 


USING 


SUB Select 


COM /Idp/ [dp 
EEP 
PRINTER IS 1 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


USING 
US ING 


USING ' 


USING 
US ING 
USING 


USING ' 


US ING 


USING ' 


USING 


USING ' 


USING 


USING ' 


Ido 


ТО Nrun 
ENTER ФЕ / је: Уа, Џо 
THEN OUTPUT AFilet:Vu.Uo 


= с (D 60 ~ оу t CO (6 — © 


"10X.""NOTE: "*,2Z2."* X-Y Pairs are stored in file ац 


ect option: ` 
Takina data or re-processing old data 
Plotting data on Loqgq-Loq '"" 

Plotting data on Linear ее 

Make cross-plot coefft file'"" 
Re-circulate uyater""" 

Purge"'"" 

TLsUÜPOD correction 5» 

Print Uo Бех 

Modify X-Y file""*"* 

Моме""" й 

Comb""" 


Fixup""" 


IF Ide-0 THEN CALL Main 
IF Idp=1 THEN CALL Plot 
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7409 
7411 
7414 
7417 
7420 
7423 
7426 
7429 
7432 
7435 
7438 
7441 
7444 
7447 
7450 
7453 
7456 
7459 
7462 
7465 
7468 
7471 
7474 
7477 
7480 
7483 
7486 
7489 
7492 
7495 
7498 
7501 
7504 
7507 
7510 
7513 
7515 
7519 
7522 
7525 
7528 
7531 
7534! 
7537! 
7540 


66),M(66 


7543 
7546 
7549 
0552 
7885 
7558 
7561 
7564 
(T) 

7567 
7570 
2523 
7576 
2579 


IF Idps2 THEN CALL Plin 

IF Idp=3 THEN CALL Coet 

IF Idp=4 THEN CALL Main 

IF Idp=S THEN CALL Purg 

IF Idp=6 THEN CALL Tdecn 

IF Idp=7 THEN CALL Uoprt 

IF Idp=8 THEN CALL Xymod 

IF Idp-9 THEN CALL Move 

IF Idp-10 THEN CALL Comb 

IF Idp=11 THEN CALL Fixup 

SUBEND 

SUB Xymod 

PRINTER IS 1 

BEEP 

INPUT “ENTER FILE NAME", Files 
ASSIGN @Filel TQ Files 

BEEP 

P "ENTER NUMBER QF X-Y PAIRS",Np 
INPUT "ENTER NEW FILE NAME",File2s 
CREATE BDAT File2s,5 

ASSIGN @File2 TO File2s 

BEEP 

INPUT "ENTER NUMBER QF X-Y PAIRS TO BE DELETED"*,Ndel 
IF Ndel=0 THEN 7492 

FOR I=1 TO Ndel 

ВЕЕР 

INPUT "ENTER DATA SET NUMBER ТО ВЕ DELETED® ,Nd(I) 
NEXT I 

РОК Js! TO Np 

ENTER @Filel:XxX,Y 

FOR I=1 TO Ndel 

IF Nd(I)=J THEN 7516 ` 
NEXT I = 
OUTPUT Ө8Ё:1е2:Х,Ү 

BRINT J.X.Y 

NEXT J 

PRINTER IS 701 

ASSIGN @Filel ТО м 

ASSIGN 3File2 TO + 

SUBEND 

SUB Move 

FILE NAME: MOVE 


E 8op(66),8(66),8(66),C(66),D(66),E(66),F(66),G(660,H(66),J(66),K(66)0,L( 


DIM Told$S(66)(14],N(66) ,Vr(66) ,Ir(66) 
BEEP 


INPUT "OLD FILE TO MOVE".D2 files 

ASSIGN @File2 TO D2 files 

ENTER 3File2;Nrun.DatesS.Ldtclt,Ldtc2,Itt 

FQR Is! TO Nrun 

ENTER ЗЕ: 1е2: Вор‹І) . То14%(1) 

ere shi lez oni 1), 8(1) CCl) .DC1).E(1) .F C1) .G01) ,H¢ 1), 501) КС) „2 01) МСП ДМ 


ENTER 3File2:Vr( D ,Ir(D 

NEXT I 

ASSIGN sFile2 TQ * 

BEEP 

INPUT “SHIFT DISK AND HIT CONTINUE". Ok 
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7542  BEEP 

7585 INPUT “INPUT BDAT SIZE",Size 

‘7588 CREATE BDAT D2_fileS.Size 

7591 ASSIGN sFilet TO D2_files 

7534 OUTPUT &Filei;Nrun,Dates.Ldtcl,Ldtc2.Itt 

7597 FÜR Tei TQ Nrun 

7600 OUTPUT ФЕ: еї: Вор 1) , То14%(1) 

76503 QUTPUT ®Filel:ACI) .B(1L),C(I) DCI) ECL). F(1),G(1) Atl) 22 Т) KOT) OD Ша 


7606 OUTPUT 8File!;Vr(D),Ir(I) 
7609 NEXT I 

7612 ASSIGN @Filel TO + 

7615! RENAME “TEST” TO D2_files 
7618 BEEP 2000..2 

7621 ВЕЕР 4000..2 

7624 ВЕЕР 4000,.2 

7627 PRINT “DATA FILE MOVED" 
7630 SUBEND | | 
7633 SUB Como 

7636! FILE NAME: COMB 


7642 DIM Emf(12) 

7645 BEEP 

7648 INPUT “OLD FILE TO FIXUPSCD2 Files 
7651 ASSIGN @File2 TO D2_files 

7654 D1! fileS-"TEST" 

7657 CREATE BDAT D1 files.30 

7660 ASSIGN *ilel TO Dt_files 

7663 ENTER ЗЕ. le2:Nrun. Dates Late! .LdadteZ itt 
7666 IF K=) THEN OUTPUT S8FilelsNrun,DateS.Ldtcl,Ldtc2,Itt 
7669 ҒОК Ізі Т0 Nrun 

7672 ENTER @File2:So0p.ToldS.Emr(*) ,МУс.Їг 
7675 OUTPUT @Filel:800,ToldS,Emt(*),Vr.Ir 
7678 NEXT I ` 
7681 А5510М 8Е:1е2 ТО » Ë 

7684! RENAME "TEST" TO D2 files 

7687 BEEP 2000..2 

7690 BEEP 4000,.2 

7693 BEEP 4000,.2 

7696 BEEP 

7689 INPUT “WANT TO ADD ANOTHER FILE (12Y,0=N)?" ,Oka 
7702 IF Oka=i THEN 

7705 Ка! 

7708  BEEP 

7711  INPUT "GIVE NEM FILE NAME",Nfiles$ 

7714 ASSIGN #File2 TO Nfiles 

7717 GOTO 7663 

¿720 ENDIF 

7723 ASSIGN ıle? TO + 

7726 SUBEND 

7729 SUB Ғіхчио 

7732! FILE: FIXUP 

7735! DATE: February 18, 1986 





7738! 
7741 DIM Emt(12) 
7744  BEEP 


7287 INPUTS “GED FIDE TU FIXUBSCD2 St Iles 
7750 ASSIGN aFile2 TO 02 files 

2055. Орје ТЕЗІ” 

7756 CREATE BDAT D1_files,20 

7753 ASSIGN @Filei TO Di_files 
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2752 
7765 
7768 
7771 
2774 
2777 
2780 
7783 
7786 
7789 
77921! 
2025 
7798 
7801 
7804 


ENTER *Fi11e2:Nrun.DateS.Ldtcl.Ldtc2.Itt 
Nrunos27 

OUTPUT GFilei:Nruno.DatesS.Ldtct.Ldtc2,.Itt 
ЁОВ Та1 10 Мгип 

ENTER 3F11e2:8op, ToldS,Emt(*),.Vr,Ir 

IF [#1 THEN 7783 

OUTPUT Е |е! ; Зор, То145. Ет? (%+) , Мг, Іг 
МЕХТ І 

ASSIGN 3File2 TO + 

ASSIGN $File! TO * 

RENAME "TEST" TO D2 files 

BEEP 2000..2 

BEEP 4000..2 

BEEP 4000,.2 

SUBEND 
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APPENDIX C 
EXAMPLES OF REPRESENTATIVE DATA RUNS 


The two following data runs (TXE222 and  GTB185) are 
representative samples of the data taken over the course of 
the investigation. The first data run for the Thermoexcel-E 
surface was for a decreasing heat-flux condition and an oil 
concentration of three percent. The second data run was for 
the  GEWA-T surface was also for a decreasing heat-flux 


condition but with an oil concentration of one percent. 


Month, date and time :05:06:11:41:20 


NOTE: Prooram name : DRP4 
isk numoer = 13 
014 file name: TXEZ222 
This data set taken on : 02:18:13:42:58 
Tube Numoer: 10 


Data Set Number ° г зок О 4 =. 3.9 02:18:13:48: 44 
ТС №: 1 2 3 4 5 5 7 8 
Temo : 8.71 9.45 3.00 8.48 9.70 8.12 10.21 8.59 
Тша Tliad Tli adé Tvaor Psat Tsumo 
8.74 2.23 2. 2.06 -1.27 -18.9 
Thetab Htube M 
6.541 1.455Е+04 9. S I9E«0A Б 
Data Set bres us - 2 Bulk Oil 2 * 3.0 02:18213:49:57 
ТС Мо: 3 4 5 5 7 8 
Temo : 8. 7 9.01 8.42 8.00 3.88 3.29 1100.33 8.51 
Twa Tliad Tliad2 Tvaor Psat Тзито 


8.55 2.24 2.22 2.09 -1.73 -18.9 
Thetao Htube 


Qdo 
6.319 1.508Е-04 9. S31E+04 


Data Set Numoer • 3. 8ulk OIl 7; = 3.9 02:18::3:50:51 

TC No: 1 2 3 4 5 6 7 8 

Temo : 8.68 8.88 8.33 7.89 9.34 8.07 10.22 8.64 
Twa ТЇїаа Tliad2 Tvapr Psat Тзито 


8.53 2.25 2:19 2.07 -t.74 -18.8 
Thetao Htuoe дар 
6.308 1.511Е+04 93.526Е+04 


Data Set Number • 4 Buik 0:11 Z ” 3.0 02:18:14:02:37 

TC No: 1 2 3 4 S 3 7 8 

Temo : 6.93 7.05 6.71 6.48 1.29 5.24 7.63 5.56 
Тша Tliad Tliad2 Tvapr Psat Tsumo 


5.56 2.25 2. e 2.14 -1.73 -18.5 
Thetao  Htubpe 
4.433 1.414Е+04 | 269E «04 


Data Set Number >° S Bulk Qil7 - 3.0 023182 4203213 

ТС Мо: 1 2 3 4 S 5 7 8 

Temo : 6.32 ТОН 6.30 6.55 7.35 6.25 7.65 6.52 
Tua Tliad Tliad2 Tvapr Psat Tsumo 


6.71 2529 2.25 2.18 -1.69 -18.4 
ТҺеїаб Нечое 


Qdo 
4.444 1.408Е+04 56. 258Е-04 


Data Set Numoer “ 6 Зий 01:17» 3.0 02:18:14:13:56 
TC No: 1 2 4 S 5 7 8 
Temo : 5.50 5.55 5.41 5.26 5.58 4.90 5.84 5.06 
Tua Tliad Tliad2 Tvapr Psat Tsumo 


5.28 2.32 2:28 2.35 -1.66 -15.3 
Thetao Htube Qdp 
2.981 1.304Е+04 3. 388E «04 
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Qata Set Number * 7 Bulk Oil £ > 3.0 
ТС Мо: ! 2 В + 5 
Temo : 5.45 5,50 5.27 5.31 5.51 
Тша Пад Tliad2 Tvapr Psat Tsumo 
5.26 2529 2.25 2.35 -1.58 -16.2 
Thetab Htube Qdo 
2.987 1.233Е-04 3. 881Е-04 


Data Set Numoer = 8 Bulk 011 Z = 3.0 
] 


TC No: 2 З 4 5 

Temo : 4.25 4.33 4.22 4.15 4.49 
Тша Tliad Tliad2 Туаөс Psat Tsump 
4.20 2.24 2.16 2.64 -1.77 -15.3 


Thetab ү: Qdo 
1.998 .176Е+04 2. 2 ASE«0A 


Data Set Number = З Bulk Qil ž * 3.0 
TC No: 1 2 3 4 5 
Temp : 4.33 4.41 4.26 4.23 4.52 


Tua Tliad Tlied2 Tvapr Psat Tsump 
4.25 2.28 2523 2:71 1.71 215,8 
Thetab шин Qdo 


2.004 .163Е-04 2.344€+04 


баға Set Numoer ? 10 Bulk Oil % = 3.0 
TC No: 1 2 3 4 5 
Temp : 3.68 3.58 3.59 3.54 3.80 
Tua Tliad Tliad2 Tvapr Psat Тзцао 
3.59 2.33 2527 2.45 -1.55 -15.4 
Thetao Htube Qdp 
1,293 1.137Е-04 1.45659E+04 


Data Set Number = 11 ш ОСС == 22 
: ; 4 
3.55 3.51 3.75 


2 Twvaor Psat Tsumo 
2.34 -1.72 -16.4 
Thetab Htube 


бао 
1.285 1.142Е+04 1.467E«0A 


No = 
оог 


Qata Set Number = 12 Buik Oil Z 3.0 
ТЕ : 1 2 g 


No 4 5 
Temp : 3.06 3.07 3.02 3.00 3.12 
Tua Tliad Tliad2 Туарт Psat Tsumo 


3401 2.43 ӚЗІ 2243 =1.69' =:7,1 
Thetab неве Qdp 
.739 1,120Е+04 3.283E«03 % 


Data Set Numoer E. 13 Bulk Gil Z = 3.0 
: Ї 


3 4 5 
Temp : 3.00 3.01 2.94 2-32 3.05 
Тша Tliad Tliad2 Tvapr Psat Tsump 
2.34 2.36 2.05 3.41 -1.77 41222 
Thetab пове дар 
. 238 .136E*04 8. 386E «03 


баса Set Number = 14 Bulk Oil 2 .. 3.0 
2 


ТС Мө: 1 5 

Temp : 2.70 2.74 2. 572 E 2.80 
Twa Tliad Tliad2 Туарг Ei. Титр 
2-/1 2.41 2.03 3.73 -1.75 -17.3 


Thetad Htubde Qdo 
.488 1.078E*94 5, 267E «03 


Qata Set Number = 15 Bulk Oil Z « 3.0 


TC No: 1 3 4 5 
Temp : 2.66 2.70 2.55 2:57 2.74 
Tua Tliad Tliad2 Tvapr Psat Tsump 


2.55 2553 2.00 3.78 -1.77 -18.0 
Thetab tuha Qdp 
. 458 .145Е+04 65, 244E«03 
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02: 182 (4214229 


5 
4.88 


7 
5.82 


8 
5.02 


02:19:14:20:33 


6 
4.07 


02:48:14221:05 
6 7 


4.10 


02:18:14:25:17 
5 7 


3.51 


loud cud 


3.44 


02:18:14:33:11 
6 7 


7 
4.58 


4.65 


3.85 


7 
3.80 


8 
4.06 


8 
4.11 


8 
3.47 


8 
3.33 


8 
2.38 3.14 2.90 
02%18%14433:43 
5 7 8 
2730 3.08 2.81 


5 
2/1 


02:18:14:40:43 
7 


2. 79 


02:18*14*41: 19 
5 


7 
2.66 22! 


8 
2.51 


8 
2.54 


Month. djate and time :95:06:11:47:44 


NOTE: Program name : ОБРА 
Disk numoer- = 10 
Old file name: GTB18S 
This data set taken on : 11:25:11:35:33 


Tube Number: 3 
Data Set Numoer = t Buik Oil % © 1.0 11:25:15:46:12 
TC No: 1 2 3 4 S 6 7 8 
Temo : 8.27 8.61 7.83 7.87 7.93 7.82 9.85 7.81 
Tua Tliqd Tliad2 Туарг Psat Tsumo 


7.37 2:23 2:28 2.24 -1.73 -16.5 
Thetab Htuoe дао 
5.739 1.581Е+04 9. 073E«04 


Data Set Number = 2 Buik Qui Z = 1.0 11:25:15:46:48 

TC No: 1 2 3 4 5 6 7 

Temo : 8.30 8.65 7.86 7.30 7.38 7.85 3.88 7.94 
Twa Tliad Tliad2 Tvaor Psat Tsumo 


8.00 2.28 2. 2 2.28 -1.68 -16.5 
Thetab Htube 
5.728 1.582Е+04 as 062E «0a 


Data Set Numoer = 3 Bulk Oil 2 = 1.0 11:25:15:54:05 
TC No: 1 2 3 4 S 6 7 8 
Temo : 6.93 7.13 6.62 6.82 6.64 6.50 7.78 6.54 
Twa Tliad Iliad2 Tvaor Psat Tsumo 
6.66 2498 2. 22% 2.77 -1.75 -17.0 


Thetab Htube 
4.448 1.348Е+04 "A 99SE+04 


Data Set Numper ° 4 Bulk Qil Z = 1.0 АРАСЫ E S4: PEN 

TC No: 1 2 3 4 S 5 

Temo : 6.35 72:15 6.61 6.62 6.62 6.49 (> "m 6. к 
Twa Tliad Tliad2 Туарг Psat Tsumo 


6.66 2.22 2: oe 2.81 -1.73 -16.3 
Thetab Htube 
4.430 1.353Е-04 = 393E«04 


Data Set Number = S Buik Qil Z = 1.0 11:25:16:03:48 

TC No: 1 2 3 4 5 6 7 8 

Temp : 6.99 6.19 5.32 5.31 6.25 5.74 6.37 S.66 
Twa Tliad Tliad2 Tvaer Psat Tsump 


5,82 2.28 2:23 3.18 =1:73 -1774 
Thetao Htuoe Odo 
3.583 1.036Е-04 3. 712Е-04 


Data Set Number = 6 Bulk Oil Z = 4.0 11:25:16:04:22 

ТС Мо: 1 2 3 4 5 6 7 8 

Temo : 6.11 6.23 5.84 5.84 5.30 5,77 6.42 5.73 
Twa ТЇїаа Tliad2 Tvaor Psat Ївшто 


5.85 25277 2.26 3.22 -1.69 -17.4 
Thetab Htube Odo 
3.588 1.034Е-04 3. 712Е+04 


Data Set Number = 7 Bulk O1i Z = 1.0 11:29: 18:10:15 
TC No: 1 2 3 4 5 5 7 8 
Temo : 5.41 6.52 6.19 5.32 S.2S 27,16 5.50 5.12 
Twa Tliad Tliad2 Tvaor Psat Tsumo 

$.24 2.24 ED 3.61 91.74  -17.7 


Thetab Нїчре 
3.020 7.320Е +03 Е 211E«04 


Data Set Number = 8 Bulk 011 % = 1.0 


11:25 008: 12* 11 
TC No: 1 2 3 4 ^ 6 7 8 
Temo : 5.25 5.32 4.98 5.13 5.06 4,95 5,29 4.30 
Twa Tliad Tliad2 Туарг Psat Tsumo 


5.04 2.14 2.05 3.56 -1.88 -17.7 
Техас Htube 


Оде 
2.949 7.480Е+03 2. 206Е +04 


Data Set Numoer = 3 Во Ок је 1.0 11525 1617 23 
TC No: 1 2 3 4 5 6 7 8 
Temo : 5.05 5.14 4.81 5.05 4.34 4.81 4.95 4.78 
Tua Tliad Tliad2 Tvaer Psat Tsumo 
4.30 2125 2.22 3.08 791.73 17.7 
Thetab Htube Odo 
2,663 5.114Е-03 1.362Е+04 
Data Set Number “ 10 Bulk 011 7 = 1.0 11:25:16:18:04 
ТС Мо: 1 2 3 4 5 6 7 8 
Temp : 5.10 5.20 4.33 5.03 4.39 4.85 5.00 4.85 
Twa Tliad 111592 Tvapr Psat Tsumo 


4.95 2:39 2227 3.44 -1.67 -17.5 
Thetab Htube 


Qdo 
2.657 5.1032+03  1.356Е»04 
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Data Set Number = 11 Bulk Qil 4 = 1.0 11:25:16:24:54 
2 5 6 8 


No: 1 3 à 7 
од 4,75 4,38 4,42 4.84 4.75 4.61 4.62 4.58 
Twa Tliad Tliad2 Tvaor Psat Tsumo 


4.66 2.30 2.23 3.15 -1.66 -16.3 
Thetab Htube Qdo 
2.360 3.273Е+03 7. 725E«03 


Data Set Number * 12 си мэ - T зше зз. 
: 1 2 
T 4.76 4.38 4.42 4.34 4.74 4.61 4.60 4,58 
Twa Tliad Т1:992 Туарг Psat Tsumo 
4,66 2.30 2.23 3.13 -1.56 -16.0 
Thetab Htubde Qdp 
2.360 3.264Е%03 7. 702E+03 
Data Set Numoer = 13 на и = E гезе 
Мо: 1 
нэ 4.45 4.54 4.09 4.60 4.39 4.34 4.36 4.28 
Twa Tliad ТЇїа442 Туарг Psat Tsump 


4.37 2527 2. с 2.45 -1.69 -14.3 
Thetab Htube 
2.094 2.361Е+03 ^s ЕТ 


Data Зее Митоег = 14 н пй: = ыд и, 
ТС №: 1 2 
Pier: 4.46 4.55 4.12 4.61 4.43 4.39 4.39 4.32 
“Twa Tliad Tliad2 Tvapr Psat Tsumo 
4.39 2.30 2223 2,44 -1.56 -14.3 


Thetab Htube дар 
2-096 2.352Е>03 4.931Е+03 


Data Set Number * 15 Bulk Qıl Z= 1.0 11:25*15:37:53 
TC No: 1 2 3 4 5 6 7 8 
Temo : 4.19 4.29 3.34 4.39 4.07 4.05 4.00 3.99 
Twa Tliad Tliad2 Tvaor Psat Tsumo 
4.09 2:23 2123 2.83 -12/3--14:8 
Thetad Htuoe Qdp 
1,863 1.582Е-03 2.946Е-03 
Data Set Number = 15 Buik 011 2 1.0 11:25:16:38:28 
TC No: 1 2 3 4 5 6 7 8 
Temp : 4.20 4:30 3.83 4.36 4.07 4.05 3.98 3.99 
Twa Tliqd Tliad2 Tvaor Psat Tsump 


4.09 2123 2. ET 2.82 -1.74 -14.8 
Thetab Htube 


1,462 1.582Е+03 >. (3 46€ +03 


* 


Data Set Number = 17 Buik Oil Z > 1.0 11:25:16:44:56 
TC No: ZI 2 3 4 5 6 7 3 
Temp : 3:97 4.04 3:67 4.17 2:77 3.31 3.57 3.74 


Twa Tliad Tliad2 Tvapr Psat Tsumo 
3.35 2.19 2419 2.71 -1.78 -15.2 
Тһасао Нечбе Qdp 

1,660 1.146Е+03 1.901Ё»03 


Data Set Number = 18 Buik Qıl Z 1.02 19225316; 45233 
TC No: 1 2 3 4 S 6 7 8 
Temp : 3:97 4.04 3.68 4.19 3.18 3.82 3.66 3.12 
Twa Tliad Tliad2 Tvaor Psat Таџто 
3.85 2413 2.19 2.71 -1.77 -15.2 
Thetab Htubde дар 
1.660 1.147E*03 1.904Е+03 


Data Set Number = 19 Bulk Qil % = 1.0 11:25:16:51: iue 
TC No: 1 2 3 å S 6 7 
Temp : 3.69 3573 3.S0 3.77 3.40 3.46 3.25 3532 
Tea Tliad Tliad2 Tvapr Psat Tsumo 
3251 2319 2:18 2.58 -1.78 -15.4 
Thetab Htube дар 
1.326 8.573Е+02  1.137E*03 
Data Set Number = 20 Buik Qil Z= 1.0 11:25:16:53:16 
TC No: 1 2 3 à S 6 7 8 
Temp : 2577 3.81 3.57 3.84 3.46 3.54 3.32 3.38 
Twa Tliad Tliad2 Tvapr Psat Tsump 


3458 2025 2.25 2.62 -1.71 -15.7 
Thetab Htube дар 
1,335 8.620Е-02 1.152Е+03 
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UNCERTAINTY ANALYSIS 


The uncertainty of the heat-transfer coefficient at ‘GE 
kW/m? апа 8 kW/m? for data runs of the GEWA-T tube (GTB185 
with one percent oil concentration) and of the Thermoexcel- 
E СГХЕ222 with three percent oil concentration) are 
analyzed below. The method of analysis is based on Kline- 
McClintock (Ref. 201] method of uncertainty analysis. 


The heat-transfer coefficient is: 


Б з (ас № Tu = т (Dene 
and 
Tuo ~= Taar = Tur - (да 1m (0,.-/0,11/С2тт Lo ЕТ (Dy Zz 
where 
h = heat-transfer coefficient 
qc = heat flux corrected for end loss 
~~ = average outer tube wall temperature 
Тват = Saturation temperature 
Та: = average inner tube wall temperature 
Qc = heat input corrected for end losses 
D, = tube internal diameter 
D2 = tube external diameter 
k = tube wall thermal conductivity 
L = length of tube's heated surface 
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To begin the analysis, a dummy variable is assigned to the 
conduction term as follows: 
DELT = иде «mweD./D,ow/( 2 T k L) ср ган 
In accordance with Kline and McClintock, the 


uncertainty of the heat-transfer coefficient is: 


Oh = 0 ас. E uw QNI. d «ні (D. 4) 
h qc Тио ~ Tsar 
_ODELT d __„ОТоњт __ 5 
Tuo = Tsar Tuo E Tsar 


The uncertainty for the conduction term is estimated since 
the logarithmic error is negligible when compared to the 
other error of the analysis. The uncertainty of the 


conduction term is as follows: 


ÓDELT 2 1-4) £ а) poi C 2 (0-5) 
DELT - Qc k L 


where Qc => Је т О з L (D6) 


The uncertainty of this term is as follows: 


` 172 
D - ба; |? + бр, \2 = ôL \? (D. 7) 
Qc Qc D> E 
Table 4 contains a listing of the computed terms of 
equations D.1 through D.7 for selected data runs for the 
GEWA-T and the Thermoexcel-E tubes at specified heat-flux 


setting of 60 kW/m? and 8 kW/m?. The algorithm in the data- 
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TABLE 4 
UNCERTAINTY ANALYSIS 


File Name GTB185 TXE222 

Heat Flux 60 kW/m2 8 kW/m2 60 kW/m2 8 kW/m2 
01175... р» 1C 9 E 
óQc. 

Qc 0.015 0.015 0.015 0.015 
ôk 

k 0.03 0.03 0.03 0.03 
óL 

G О 8025 0.0025 0.0025 0.0025 
DELT 0.486 0.064 |. 0.486 0.064 
óDELT 

DELT 0.069 0.520 0.069 0.526 
TOM OEC 6.845 4.678 6.904 2.964 
Tu ТСС) 0.49 0.065 0.49 0.065 
Tete. (CC) 22215 2.295 2.265 2.205 
Tones AUR) O. 1 5 O. 1 O. 1 
оаа 
Е 0.111 0.028 0311! 0.088 
£ oDELT 
Тото O. 008 0.014 0.008 0.046 
са Таж _ 
ТА Тел" 0.023 0.042 0.023 0.135 
бас. 

ас 0.02 0.02 0.02 0.02 
oh 

h 0.115 0.056 0.115 0.169 
h(W/m2 K) 13480 3273 14140 11200 
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reduction program used this approach and assumed that all 
the eight wall thermocouples were used to calculate the 


heat-transfer coefficient. 
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